A0-A174  691  HOMOGENIZING  SURFACE  AND  SATELLITE  OBSERVATIONS  OF  '  1/1 

CLOUDCU)  EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH  ANO 
DEVELOPMENT  FPO  NEU  A  HENDERSON-SELLERS  ET  AL 
UNCLASSIFIED  28  NOV  86  EOARD-TR-87-81  AFOSR-85-0299  F/G  4/1  NL 


MIC  FlirCOPY  AD-A174  691 


UJHMPJ  PJWJ  VU 


Contract  /  Grant  Number  AF0SR-&5-0299 


HOMOGENIZING  SURFACE  AND  SATELLITE  OBSERVATIONS  OF  CLOUD 


A.  Henderson-Sellers  and  A.  Goodman 
Department  of  Geography, 

University  of  Liverpool, 

Liverpool,  U.k. 


28  November  1 9&6 


Interim  Scientific  Report  ,  1  October  19&5 

5 0  September  198b 


Approved  for  public  release;  distribution  unlimited 


Prepared  for 


EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH  AND  DEVELOPMENT 
London,  England. 


^  DTIC 

FLECTEf* 
DEC  1  1986  I 


%  /£  01  OJ 


8  6  12  01  015 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


% 


-(WRi*  ■  '* J 


-7  3  V 


This  report  has  been  reviewed  by  the  EOARD  Information  Office  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

CW  ft. 

OWEN  R.  COTE’ 

Chief,  Geophysics  and  Space 

- 

ROBERT  C.  WINN,  Lt  Colonel,  USAF 
Chief  Scientist 


»  ,4  f  •  r 

iV  c/ 


[l* 


REPORT  DOCUMENTATION  PAGE 


I.  Report  Number 


4.  Title  (and  Subtitle) 

Homogenizing  surface  and  satellite 
observations  of  cloud. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  Recipient  s  Catalog  Nunber 


5.  Type  of  Report  &  Period  Covered 

Interim,  1  October  85  - 
30  September  86. 


6.  Performing  Org.  Report  Number 


7.  Author (s)  8.  Contract  or  Grant  Number 

Ann  Henderson-Sellers  and  Alan  Goodman  AFSOR-85-0299 


9.  Performing  Organization  Name  and  Address 

Department  of  Geography, 

University  of  Liverpool, 

P.0.  Box  147, 

Liverpool,  L69.  3BX. 


11.  Controlling  Office  Name  and  Address 
EOARD 

223/231  Old  Marylebone  Road, 
London,  NV1 ,  5TH  U.K. 


14.  Monitoring  Agency  Name  and  Address 


10.  Program  Element,  Project,  Task 
Area  &  Work  Unit  Numbers 
Program  element:  6ll02F 
Project/Task:  2301/Dl 


12.  Report  Date 
28  November  1986 


13.  Number  of  Pages 
82 


15. 


16.  &  17.  Distribution  Statement 

Approved  for  public  release;  distribution  unlimited. 


18.  Supplementary  Notes 


19.  Key  Words 

Cloud  retrieval,  satellite  radiances,  surface  observations  of 
clouds . 


20.  Abstract 

^  The  major  types  of  cloud  retrieval  algorithms  are  reviewed  with 
special  emphasis  being  placed  upon  recent  { i  re-» — post  1 98-1/82)  Q 
developments  and  novel  techniques.  Satellite/based  retrieval 
algorithms  can  be  grouped  into  three  classes:  threshold  methods, 
statistical  procedures  and  radiative  transfer  techniques  although 
each  algorithm  depends  upon  implicit,  if  not  overt,  inversion  of 
the  radiative  transfer  equation.  The  fourth  type  of  retrieval 
differs  very  considerably  from  the  satellite-based  techniques  as 
it  depends  upon  surface-based  measurements  which  are  generally, 
but  not  always,  human  rather  than  computer  based.  There  is  a 
tendency  to  assume  that  surface-based  observations,  especially 
of  total  cloud,  are  ^correct'*  but  that  they  may  differ  from 
satellite*based  retrievals  because  of  the  differences  in  viewing  — 
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geometry.  Actually  the  literature  reveals  surprisingly  few 
intercomparison  studies.  None  of  the  satellite-based  techniques 
have  yet  been  well  validated  in  a  variety  of  situations  while 
surface-based  observations  are  made  in  all  terrain  and  climate 
regimes.  Overall  there  seems  to  be  good  reason  to  believe  that 
surfacef-based  observations  of  cloud  amount  and  especially  of  low 
cloud  amount,  character  and  base  height  can  add  significantly 
to  satellite-rbased  retrieval  algorithms.  ’}i.  1  ' 
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1_  Introduction 

This  report  comprises  the  first  stage  of  a  Ph.D.  study.  The  aim 
is  to  review  the  basis  of  cloud  observation  techniques  laying 
particular  emphasis  on  satellite-based  retrievals.  Since  the  early 
19&0s  there  have  been  a  number  of  developments  in  remote  sensing  of 
cloua  characteristics  and  these  are  reviewed  here.  It  is  clear  that 
whilst  many  of  the  more  novel  techniques  offer  additional  insight  into 
cloud  features  none  of  them  have  yet  been  adequately  validated. 

The  li.S.  Air  Force's  nephana  lyses,  3b  Neph  and  RT  Neph  have  been 
excluded  intentionally  from  this  review  for  two  reasons:  (i)  they  are 
based  on  a  range  of  data  from  a  variety  of  sources  and  (ii)  an 
adequate  description  of  the  newer,  current  algorithm  (RT  Neph)  is 
still  awaited.  Despite  this  exclusion  our  conclusions  (see  Sections  3 
and  4)  suggest  that  RT  Neph  may  well  be  benefitting  greatly  from  the 
inclusion  of  conventional  as  well  as  satellite-baseo  observations. 

With  the  advent  of  satellite  technology,  we  have  been  provided 
with  the  facility  to  observe  clouds  from  above  as  well  as  below.  Lue 
mainly  to  their  vastly  superior  areal  coverage  the  last  two  decaoes 
have  witnessed  a  sustained  effort  in  utilising  remotely  sensed 
radiation  from  various  regions  of  the  electromagnetic  spectrum  to 
derive  cloud  parameters  relevant  to  meteorological  and  climatological 
studies.  The  human  eye,  though  still  employed  in  a  dal  1  y 
observational  sense  has  been  forced  to  take  a  secondary  role  compared 
to  that  of  the  satellites.  Ihe  development  and  use  of  sopnisticated 
climate  models  has  taken  into  consideration  the  radiative  effects  of 
the  global  cloud  field  ano,  along  with  increasing  awareness  of  the 
significant  role  the  various  cloud  types  play  in  determining  the 


Earth's  radiation  budget,  has  contributed  to  the  neea  for  global  clouo 
coverage  information  which  can,  of  course,  only  be  achieved  from 
spaceborne  sensors,  there  being  insufficient  resources  available  to 
deploy  the  required  density  of  surface  observers,  especially  over  the 
oceans# 

The  International  Satellite  Cloud  Climatology  Project  (1SCCP) 
(korld  Climate  Programme  (WCP),  19 k>2;  Schiffer  and  Rossow,  19&3,  19&b) 
aims  to  provide  users  with  a  uniform  global  radiance  climatology  from 
which  specific  cloud  parameters  will  be  extracted.  The  method  of 
extraction  (i.e.  choice  of  particular  cloud  retrieval  algorithm)  is  of 
crucial  importance.  As  will  become  evident,  there  is  at  present  no 
single  retrieval  algorithm  capable  of  performing  all  the  tasks 
required  of  it  and,  equally  important,  a  lack  of  'truth'  data  sets 
against  which  results  can  be  compared  (Rossow  e_t  a  1# ,  19bb).  For 
ISCCP  to  be  of  sustained  beneficial  use  the  cloud  retrieval  (as  well 
as  being  as  accurate  as  possible),  will  need  to  be  very  precise  (Table 
1  )•  In  this  case  the  accuracy  of  retrieval  refers  to  the  minimising 
of  systematic  errors  in  the  reception  and  calibration  of  the  satellite 
data#  The  high  precision  requirement  defines  the  confidence  limits  of 
the  estimated  clouo  parameters.  Many  recent  algorithm  developments 
have  taken  place,  in  many  cases  building  on  previous  work  to  increase 
the  sensitivity  towards  detection  of  specific  cloua  systems.  In  the 
following  sections  present-uay  knowledge  of  satellite  cloud  retrieval 
is  reviewed,  the  main  problems  analysed  and  the  feasibility  of  merging 
surface  observations  into  the  retrieved  cloud  fields  examined  in  the 
light  of  recent  work  (e.g.  Sfeze  e_t  ah,  1 9 b 6 )* 

khen  the  first  sate  1 1  i.te-oeri.  ved  images  were  made  available  to 


Table  J_  Data  specification  for  the  International 
Cl imatology  Project  (after  Rossow  e_t  a  1.,  19b5) 


Parameters  -  Spatial  and  temporal  averages  and  va 
statistical  measure  of  the  shape  of  the  temporal 
required  for  each  of  the  following  parameters. 


Amounts 

g 

Total  cloud  amount  (fraction) 

Cirrus  cloud  amount  (fraction) 

Middle  cloud  amount  (fraction) 

Low  cloud  amount  (fraction) 

Deep  convective  cloud  amount  (fraction) 

Height 

Cirrus  cloud-top  height  (km) 

Middle  level  cloud-top  height  (km) 

Low-level  cloud-top  height  (km) 

Deep  convective  cloud-top  height  (km) 

* 

Cloua-lop  Temperature  ( k)  for  each  cloud  category 
Cloud  Optical  Depth 
Cloud  Size  Distribution 

Average  Narrow  Band  Radiances  (VIS  and  IR) 

Spatial  Averaging  -  The  information  is  to  be 
approximately  250  km  by  250  km  boxes 

Time  Samp  1  ing  -  Every  3  hours  i.e.  b  times  a  day, 
synoptic  observation  times 

Time  Averaging  -  The  global  cloud  climatology  wil 
averages  for  each  of  the  6  observing  times  per 

Length  of  Time  Series  -  5  years 
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+0.03 
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+0.05 
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+  1.00 
+0.50 
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+  1.00 
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day 
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users,  meteorologists  were  faced  with  a  choice  of  either  carrying  out 
a  manual  analysis  of  the  photographic  image  (Clapp,  1964''  or  devising 
an  automated  technique  to  process  the  digitized  data,  in  order  to 
ascertain,  for  example,  cloud  amount  ( A r k i n g ,  1964).  Manual 
nephanalyses  were  found  to  consume  many  man-hours  of  time  whilst  at 
the  same  time  incorporating  an  inevitable  degree  of  subjectivity. 
Automated  processing,  although  unsuccessful  initially  (Gooshall,  1970) 
is  now,  with  the  aid  of  modern  computers,  a  much  more  efficient  means 
of  analysis,  especially  for  quantitative  determination.  ISCCP  will 
eventually  provide  over  1600  days  worth  of  global  data:  a  data  archive 
clearly  beyond  the  ability  of  manual  analysts. 

2  Cloud  Retrieval  Algorithms 


All  cloud  retrieval  algorithms  comprise  two  basic  steps,  namely 
initial  detection  of  cloud  followed  by  quantitative  analysis  of  that 
cloud.  The  first  step  involves  separating  observed  radiances  into 
either  a  clear  or  cloudy  category.  The  methodology  employed  in  the 
separation  defines  a  particular  type  of  algorithm.  The  analysis  step 
concerns  the  determination  of  cloud  properties  from  the  observed 


radiances  and  may  vary  from  the  fitting  of  complex  radiative  transfer 
models  in  order  to  obtain  parameterised  cloud  properties  to  the 


straightforward  task  of  counting  clouoy  pixels  to  determine  cloud  cover 
fraction.  Implicit  in  any  algorithm  is  a  set  of  assumptions  (in  effect 
a  radiative  model)  characterising  the  relationship  between  the  Earth's 
surface,  the  atmosphere  and  clouds  with  the  satellite  measurements. 


Algorithms  may  be  distinguished  by  the  type  and  complexity  of  the 
inherent  radiative  model  but  all  are  designed  to  retrieve  cloud 
properties,  the  retrieval  being  dependent  upon  the  methodology  used. 
The  accepted  means  of  di. sti. notion  is,  however,  the  application  of 


detection  and  analysis  to  either  groups  of  image  pixels  or  individual 
image  pixels.  This  latter  type  of  algorithm  is  termed  a  threshold 
method  whilst  the  former  type  is  known  as  a  statistical  method.  In 
addition  a  third  type  of  algorithm,  the  radiative  transfer  technique 
can  be  arbitrarily  defined  by  virtue  of  its  differing  approach  to 
retrieval.  Each  category  ana  its  applications  are  reviewed  in  turn 
and,  where  possible,  illustrated  with  reference  to  NOAA-7  polar 
orbiting  data  for  Western  Europe  on  10th  August  1  965.  The  data 
comprise  three  arrays  of  51k  x  51k  image  pixels  covering  the  latitude 
belt  42°  to  57°  North.  The  major  weather  systems  present  are  a 
depression  in  the  Eay  of  Biscay  with  associated  convective  disturbances 
on  its  eastern  flank  and  an  anticyclone  with  clearer  skies  over 
Southern  Scandinavia.  The  images  (Figure  1(a)  is  channel  1  [0.56  -0.6b 
j*m]  and  Figure  1(b)  is  channel  5  [11.5  -  Ik. 5  Jim]),  which  are  for  1425 
GhT  on  the  10th  August  1963,  and  the  surface  synoptic  chart  for  midday 
on  the  10th  August  (Figure  1(c))  indicate  a  range  of  features:  fog  in 
the  North  Sea,  multi-level  cloud  over  France  and  some  clear  land  and 
ocean  areas. 

2*1  Threshold  Techniques  for  Cloud  Retrieval 

The  threshold  method  treats  each  image  pixel  separately, 
assigning  to  each  pixel  or  field  of  view  a  completely  clear  or  cloudy 
label  according  to  the  magnitude  of  the  observed  radiance  relative  to 
the  predetermined  threshold  level  i.e.  the  pixel  cloud  fraction, 
(referred  to  hereafter  as  f)  equals  either  0  or  1.  This  definition 
(Arking  and  Childs,  1965)  appears  to  contradict  that  given  in  the 
1SCCP  preliminary  implementation  plan  (WCP,  19621  which  states  that 
the  cloud  fraction  i.s  linearly  related  to  the  difference  between  the 
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observed  ana  threshold  radiance,  implying  values  of  f  intermediate 
between  0  and  1.  however,  Godshall  (1570),  reporting  on  the  previous 
work  of  Barnes  ano  Chang  (1966)  implied  a  definition  consistent  with 
the  former  (i.e.  f  =  0  or  1)  ano  applications  reported  since  1962 
(Liljas  1564,  Wielicki  ano  Welch,  1566)  confirm  this  view* 

A  single  channel  threshold  was  originally  used  (e.g.  Saunders 
and  hunt,  1 9  &  0  )  which  involved  no  auxiliary  data  to  determine  cloud 
cover,  this  method  is  termed  the  ’gross  clouo  check'  in  Saunoers 
(1966).  however,  bispectral  or  double-channel  thresholds  are  now  also 
documented  (Minnis  and  harrison,  1984a, b,c  ano  see  also  Ackerman  ano 
Cox,  15&1).  The  ISCCP  pilot  study  of  1561  (World  Climate  Programme 
1962),  aiming  to  establish  the  relative  effectiveness  of  the  then 
available  algorithms  included  both  single  channel  ano  bispectral 
thresholds  in  its  evaluation.  The  single  channel  version  defines  the 
threshold  as  some  increment  in  either  visible  reflectance  ( AR)  or 
infrared  brightness  temperature  (AT)  from  a  clear  sky  value  (hs,Ts) 
so  that  for  visible  (infrared)  radiances,  a  pixel  is  labelled  cloudy 
if  R>hs  +  Aft  (or  T<TS* AT).  The  cloud  cover  fraction  can  then  be 
obtained  by  counting  up  all  the  'cloudy'  pixels  whilst  other 
parameters  e.g.  visible  albedo,  cloud  top  temperature  can  be  derived 
from  clouo  radiative  models  and  the  cloudy  radiance  values.  The  clear 
sky  and  cloudy  radiances  form  a  monotonic  distribution,  the  clear  sky 
radiance  corresponding  to  minimum  reflectance  (visible  channel)  or 
maximum  brightness  temperature  (infrared  channel).  Time  records  of 
the  extreme  values  for  each  pixel  are  examined  and  1  b  or  30  day 
composite  images  can  be  constructed  (e.g.  see  Seze  ano  Lesbois,  1966) 
where  a  pixel's  minimum  reflectance,  say,  is  assumed  to  represent  the 
clear  sky  reflectance  for  that  period.  In  certain  situations  minimum 


visible  reflectance  aoes  not  coincide  temporally  kith  maximum 
brightness  temperature  (S4ze  and  Desbois,  1 9 & 6 ) •  Ihe  time  period 
should  be  sufficiently  long  so  as  to  account  for  any  variation  of 
surface  properties  and  in  the  case  of  geostationary  satellites  it  can 
compensate  for  the  effects  of  viewing  geometry  (which  do  not  alter  in 
time).  Spatial  filtering  of  the  composite  image  can  remove  any  noise 
due  to  cloua  contamination  ana  then  the  threshola  increments, 
typically  3%  reflectance  and  6  R.  brightness  temperature  can  be 
applied# 

Alternative  methods  of  setting  the  threshold  exist:  visual 
inspection  of  the  one-dimensional  frequency  histogram  of  either 
visible  reflectances  (Wielicki  ana  Welch,  19&6;  Saunaers,  19&6)  or 
spatial  variations  in  the  raaiance  field.  In  the  former  case,  appliea 
by  Wielicki  ana  Welch  ( 1 9 b 6 )  to  fair  weather  cumulus  fielas,  the 
reflectance  histogram  (Figure  2)  displayea  a  highly  uniform 
distribution  of  cloua  reflectance  ana  the  cloua  fraction  was  strongly 
related  to  the  selectea  threshola  as  shown  in  Figure  5»  In  this  case 
it  was  concluaea  that  many  small  clouas  ( ~ 1 00m  in  size)  ha  a  been 
observed  (because  of  the  fact  that  reflectance  is  observed  to  increase 
notably  with  cloua  size)  ana  that  reflectance  across  inaiviaual  clouds 
was  highly  irregular.  Alternatively,  Saunders  ( 1 9 b 6 )  applied  a 
'aynamically  controlled'  threshold  using  several  inter-related 
parameters  to  separate  any  cloud  free  peaks  present  from  remaining 
clouay  radiances  (Figure  4).  Examination  of  the  histograms  in  Figure 
5  shows  that  in  the  visible  case  (Figure  5(a))  values  of  the  threshold 
parameters  m,  n  ana  T  coula  be  tuned  to  give  an  accurate  cloua 
fraction  estimate  whilst  in  the  infrared  example  (Figure  5 ( b > >  the 
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Histogram  of  Landsat  digital  count.  (  i  •  e  .  , 
reflectance)  for  a  western  Atlantic  cumulus  cloud  field. 
The  curve  is  plotted  for  all  127  digital  count  values  but 
for  simplicity  selected  values  are  shown.  The  equivalent 
nadir  lambertian  reflectance  scale  is  given  at  the  top  of 
the  figure  (after  Wielicki  and  Welch,  1986). 


Nadir  reflectance,  percent 
2  6  10  14  U 


0  10  20  30  40 

Cloud  threshold,  digital  counts 


Cloud  fraction  as  a  function  of  cloud  reflectanc 
threshold  for  the  western  Atlantic  cloud  field  (afte 
Wielicki  and  Welch,  1  966). 


Examples  of  (A)  visible  (0.5&  -  0.68^jm)  radiance  and 
(fi)  11  micrometer  brightness  temperature'  histograms  with 
cloud  free  peaks  (after  Saunders,  1986).  The  computed 
dynamic  threshold  is  given  as  T  which  is  determined  as 
follows  from  visible  (infrared)  histograms:  first 
calculate  the  number  of  radiances  that  constitute  the 
peak  and  the  value  of  I  |c>  If  the  peak  is  considered  to 
be  significant  then  (i)  de ' ermine  the  lower  and  upper 
limits  of  the  histogram  (i^  and  I2  respectively),  (ii)  If 
I  1  ( 1 2  ^  lies  within  'm'  counts  of  the  peak  and 

^pk^m^^^min^'  then  the  peak  is  assumed  to  be  cloud 
free  with  tne  threshold  T  set  at  'n'  counts  above  I  ^ . 
Values  of  the  parameters,  m,  n,  I  m  _  and  T  •  are  giifen 
in  Saunders  (1986).  max  min 
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Figure  5(a)  One  dimensional  frequency  histogram  for  AVHRR  channel  1 
(0.56  -  0.68  uni)  taken  from  a  sub-scene,  frame  B  (47° 
-52°)  of  Figure  1(a). 
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close  proximity  of  the  cloua  free  and  cloudy  peaks  means  that  small 
changes  in  m,n  ana  T  will  proauce  larger  fluctuations  in  the  derivea 
cloud  fraction.  The  latter  case  for  spatial  variations  assumes  that 
areas  larger  than  the  pixel  resolution  possess  homogeneous  surface 
properties  and  are  never  completely  cloud-coverea.  1  he  retrieval 
technique  is  therefore  limitea  by  surface  variation  over  areas  larger 
than  an  average  cloud  size  which  may  be  as  large  as  250  km  x  250  km. 
One  inherent  problem  with  this  method  is  that  it  does  not  offer  a 
method  of  identifying  apparent  cloua  variations  resulting  from 
variations  in  viewing  geometry  or  atmospheric  conditions. 

The  major  drawbacks  concerning  the  operational  use  of  threshold 
algorithms  relate  to  situations  of  partially  clouaea  pixels  and  the 
definition  of  the  clear  sky  radiance  values.  The  problem  of  small 
clouds  failing  to  fill  up  individual  fields  of  view  has  been  addressed 
(Coakley  and  bretherton,  1  9  b  2 ;  hossow  et  a 1« ,  19b5)  ana  might  well 
apply  in  the  case  of,  say,  fair  weather  cumulus  (hielicki.  ana  helch, 
19bb).  If  a  threshold  is  selected  to  counter  errors  in  mean  cloua 
cover  fraction  for  one  cloua  type,  considerable  biases  may  well  then 
arise  for  other  cloua  types  (Coakley  and  bretherton,  1 9 b 2 ) .  In 
addition,  when  only  part  of  the  field  of  view  is  covered,  cola,  thick 
clouds  will  appear  warmer  ana  thinner  whilst  warm,  thin  clouds  might, 
in  an  extreme  case,  appear  as  a  clear  view  to  the  surface  thus 
introducing  serious  errors  in  determination  of  cloua  amount,  cloud 
optical  properties  and  cloud  top  temperature.  Coakley  ana  bretherton 
( 1 9 b 2 )  apply  the  findings  of  Shenk  and  Salomonson  (1972)  to  hignlight 
the  risks  associated  with  thresholds  when  the  satellite's  field  of 
view  is  smaller  than  the  typical  cloua  size  by  as  much  as  two  orders 
of  magnitude.  Plotting  cloud  cover  errors  against  the  inverse  square 


root  of  clouo  area  to  scanning  pixel  size  yields  an  approximate  linear 
relation  showing  that  the  errors  decrease  inversely  with  the  square 
root  of  cloud  area  (and  hence  cloud  perimeter)  thus  explaining  the 
likelihood  of  large  thresho lo-baseo  estimation  errors  for  relatively 
small  cloud  size  distributions. 

The  logical  step  to  offset  errors  due  to  partial  pixel  coverage 
is  to  reduce  the  pixel  resolution  as  far  as  possible,  thus  increasing 
the  likelihood  of  completely  clear  or  clouoy  pixels.  Unfortunately 
for  global  analysis  pixel  resolution  is  stringently  controlled  by  the 
volume  of  resulting  data  so  that  whilst  the  Lanosat  thematic  mapper 
(T  M )  has  a  resolution  of  30  m  ano  the  Landsat  mu  Itispectra  1  scanner 
(NSS)  a  resolution  of  57  m,  as  opposed  to  hOAA's  advanced  very  high 
resolution  radiometer's  (AVhFiR)  value  of  1.1  km,  1SCCP  will  process 
only  geostationary  and  AVhkh  data  which,  will  in  turn,  require 
subsequent  compression.  Ultimately,  the  accuracy  of  clouo  cover 
determination  rests  with  the  success  of  the  threshold  technique  in 
partitioning  partially  cloudy  ano  clear  pixels  which  in  turn  depends 
on  the  magnitude  of  the  threshold  (assuming  clear  sky  radiances  are 
known  accurately)  ano  on  the  methodology  used  to  apply  it. 

The  determination  of  clear  sky  radiances  requires  that  at 
sometime  curing  the  time  period  examined,  cloud  free  conditions 
prevail  over  each  pixel.  Providing  that  surface  and  atmospheric 
properties  remain  constant  throughout  the  time  period,  the 
distributions  of  h  and  Tg  over  the  pixels  may  be  accurately  obtained 
and  the  clouds  become  discernible  through  a  single  spectral  channel, 
however,  suoden  or  gradual  spatial  and/or  temporal  variations  in  the 
surface  or  atmospheric  properties  such  as  a  heavy  fresh  snowfall  or  an 
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increase  in  atmospheric  turbidity  aue  to  an  influx  of  desert  aerosol 
can  render  the  clear  sky  radiances  inaccurate.  For  example  hossow  et 


g  ]•• »  ( 1 9  b  5 )  describe  such  a  sequence  of  events  that  occurred  during 
the  ISCCP  pilot  study  test  period.  Assuming  reliable  clear  sky  values 
the  use  of  small  thresholds  can  allow  for  the  detection  of  clouds 
whose  spectral  signatures  closely  resemble  that  of  a  clear  sky,  e.g. 
low  level  clouds  detected  by  use  of  a  discriminating  visible  channel 
threshold  (Saunders  1986). 

The  presence  throughout  the  world  of  spatial  inhomogeneity  and 
temporal  variation  of  the  surface,  means  that  the  use  of  fixea 
thresholds  for  analysis  of  global  radiance  data  is  completely 
impractical.  Ihresholos  would  have  to  be  'tuned'  to  the  situation  in 
hand  using  auxiliary  data,  hot  surprisingly,  one  of  the  major 
findings  of  the  ISCCP  pilot  study  was  the  urgent  need  for  improvement 
in  the  techniques  to  specify  clear  sky  radiances.  Future  improvements 
in  satellite  cloud  retrieval  are  likely  to  arise  directly  as  a  result 
of  a  superior  determination  of  the  clear  sky  scene. 

Application 

Despite  the  attendant  problems  outlined  above,  threshold-based 
cloud  retrieval  techniques  have  been  successfully  used  in  several 
instances  (Li  1 j a s  1984;  h inn is  and  harri son,  1984a, b , c ;  Saunoers 
198b).  The  accurate  determination  of  sea  surface  temperature  via 
remote  sensing  may  be  used  as  data  input  to  climate  and  forecast 
models.  Conventional  sounding  techni. ques  rely  on  cloud  free 
conditions  for  operation  which,  during  periods  of  persistent  scattered 
or  total  cover  may  be  an  unacceptably  long  time  in  forthcoming.  The 
dynamic  visible  threshold  method  (Saunders  198b)  has  been  successfully 
utilised  i.n  daytime  hours  as  a  component  of  a  larger  algorithm  known 


as  the  spatial  coherence  visible  method  to  provide  accurate  sea 
surface  temperature  values  by  comparison  with  ship  measurements,  here 
a  fixed  threshold  (the  gross  cloud  check)  immediately  rejects 
obviously  clouoy  pixels  whilst  the  dynamic  threshold,  operated  in 
combination  with  the  spatial  coherence  method  (see  Section  2.3) 
distinguishes  sea  surface  from  low  cloud  whilst  remaining  unaffected 
by  any  strong  sea  surface  temperature  gradients  in  the  vicini.  ty  of 
frontal  regions. 

In  attempting  to  estimate  diurnal  variations  of  cloudiness  in 
association  with  clear  sky  parameters,  the  'hybrid  bispectral' 
threshold  technique  (hinnis  ana  Harrison,  1  5  b  4  a  ,  b  ,  c  )  has  been 
introduced  ano  applied  to  subregions  of  GOES-East  images.  At  night 
the  method  reduces  to  a  simple  infrared  threshold  but  the  daytime 
visible  radiances  are  employed  to  tune  the  night-time  threshold  and 
thus  the  consistency  is  preserved  throughout  the  diurnal  cycle  (Minnis 
and  Harrison  19&4a).  The  clear  sky  visible  reflectances  are  obtained 
from  time  records  of  minimum  reflectance  for  observations  at  a  time 
close  to  midday,  R  then  being  assumed  constant  throughout  the  time 
period.  Empirically  derived  reflectance  models  which  take  account  of 
surface  properties,  solar  zenith  angle  and  the  satellite  viewing 
geometry  enable  determination  of  hs  at  other  times.  Average  values  of 
Rs  for  each  subregion  of  250  km  x  250  km  are  obtained  wi  th  /\K  set  at 
approximately  1%  reflectance.  Ihe  brightness  temperature  threshold  is 
then  determined  by  summing  (for  decreasing  values  of  'I'*  all  the  pixels 
for  which  R<RS  +  A  ^  until  their  average  value  is  equal  to  T,  the 
average  of  all  the  pi. xels  where  R<RS  +  A  R.  The  final  value  of  T  used 
in  the  summation  defines  T  -  Ai  and  all  pi. xels  colder  than  this 


value  are  labelled  cloudy  (irrespective  of  reflectance''.  In  this  way 
despite  the  inclusion  of  slightly  cloudy  pixels  or  those  containing 
cirrus  cloud  (poorly  detected  by  visible  thresholds)  which  may  cause  T 
<TS  the  effect  of  counting  partially  cloudy  pixels  as  totally  cloud  is 
offset.  In  the  next  stage  all  ft  ,  ft  and  I  values  are  normalised,  thus 
removing  any  angular  dependence  or  atmospheric  effects.  Empirical 
reflectance  models  for  land,  cloud,  ocean  and  a  clear  atmosphere  are 
applied  to  the  reflectances  to  normalise  them  with  respect  to  the 
solar  z e n i t h  angle,  the  satellite  viewing  zenith  angle  ano  the 
satellite  viewing  azimuth  angle,  whilst  application  of  seasonal  and 
latitude-dependent  limb  darkening  functions  helps  eliminate  any  water 
vapour  absorption  effects  from  the  T  values  by  normalising  to  the 
overhead  sate  1 1 ite- v iew ing  zenith  angle  of  zero  degrees  and  thus 
correcting  for  slanted  paths,  reducing  the  amount  of  atmospheric 
attenuation  possible  by  water  vapour  absorption,  (Einnis  ano  harrison 
The  clear  pixels  perniit  surface  temperature,  averaged  over  a 
subregion  to  be  estimated,  ano  the  cloud  optical  properties  are  deduced 
from  R,  T,  Rs,  Ts  and  f,  for  each  pixel.  These  properties  can  be 
assigned  to  a  low,  middle  or  high  clouo  category  by  defining,  for  each 
category,  a  range  of  temperatures  interpolated  from  Ts  assuming  a 
uniform  lapse  rate.  The  results  of  the  hinnis  and  harrison  ( 1 Sb^a, b,ci 
application  of  the  bispectral  threshold  to  monthly  GOES  data  could  be 
compared  with  the  monthly  clouo  distributions  obtained  from  long  term 
climatological  values  (e.g.  hahn  e_t  a  1 . ,  196*0  as  a  possible 
' validation  * . 

The  potential  uses  of  automated  clouo  classification  include 
climate  modelling,  numerical  weather  forecasting,  solar  energy 
utilisation  ano  agricultural  applications.  Ihe  European  Space 


Operations  Centre  (ELOCi  already  operates  such  a  procedure  for  mapping 
areal  distributions  of  the  cloud  top  level  whilst  the  Swedish  Space 
Corporation  (SSC)  employs  an  automated  classification  algorithm  for 
use  in  forecasting,  describee  by  Liljas  ( 1 9  b  4  > »  The  method  is  known 
as  mu  1 1 ispec t ra 1  box  classification  and  makes  use  of  four  AVhRR 
channels  (1,  2,  3  ana  0.5b  -  0.b9  jJ m,  0.73  -  1.10  jiru,  3.55  -  3.53 
^im,  10.5  -  11.5  ^im)  to  examine  the  radiative  properties  of  clouas  ana 
land,  ocean  and  cryospheric  surfaces  as  a  function  of  their  optical 
depth,  height,  density  and  thermodynamic  phase.  Li  1 j  a  s  ( 1 5  b 4  ) 
describes  the  spectral  characteristics  of  each  common  cloud  type  for 
the  channels  used  (Figure  b),  and  demonstrates  how  the  spectral 
similarities  between  snow,  mist,  sunglint  and  clouas  in  channels  1  ana 
4  can  be  overcome  by  the  use  of  channel  3  (Figure  7).  where  the 
application  of  channel  3  brightness  temperature  minus  channel  4 
brightness  temperature,  normalised  by  division  by  the  sine  of  the  solar 
elevation  angle,  E.  provides  a  measure  of  the  reflectance  of  solar 
radiation  in  channel  3  (Table  2)  (Liljas,  1564).  The  implementation  of 
fixed  thresholds  in  the  n-dimensional  radiance  intensity  space  (where  n 
represents  the  number  of  spectral  channels  usea)  i s  used  to  define 
'boxes'  for  16  cloud  and  surface  types  (Table  3)  for  the  summer  period 
Fiay  to  September.  The  method  shows  considerable  agreement  in 
comparison  with  synoptic  surface  observations,  although  a  comparison 
with  different  (i.e.  non- thresho Id)  types  of  algorithm,  whilst  not  yet 
attempted  would  prove  a  useful  exercise.  Unfortunately  Liljas  (19641 
fails  to  provide  information  concerning  the  sensitivity  of  the  applied 
thresholds  to  changes  in  season  (ana  hence  surface  properties)  although 
some  level  of  tuning  seems  likely  to  be  necessary.  however,  from  the 


0 


gure  6^ 


40 


60 


W  mam va,  u?i 


Position  of  various  surfaces  in  visible  (channel  1, 
0.5b  -  0.68  um)  and  infrared  (channel  4,  10.5  -11.5  p®) 
intensity  space  from  AVHRR  data  (after  Liljas,  1984). 
Standard  U.K.  Meteorological  Office  cloud  symbols  are 
used  . 
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Figure  7 


Positions  of  different  reflecting  surfaces  in  two- 
dimensional  intensity  space.  The  vertical  axis  is  the 
brightness  temperature  difference  (°C)  in  channels  3 
(3*55  -  3*93  jim)  and  4  (10.5  -  11.5  jie)  normalised  by 
division  by  the  sine  of  the  solar  elevation  angle,  E,  to 
an  overhead  sun  value.  The  horizontal  axis  is  the 
normalised  cloud  reflectance  in  channel  1  (0.58  -  0.68 
jim)  (after  Liljas,  1984). 


in  AVhRR  channels  3  and  4.  Mean  value  and  extremes  for  May  to 

September  of  well  defined  surfaces 

(after  Liljas,  1964) 

• 

Cloud  /  surface  type 

Albedo  (Tv 

-T a ) /sin  E  (-C) 

[ 

extreme  -  mean 

extreme  ] 

Thin  stratus/fog 

0.20-0*40-0.52 

20-27-36 

Stratocumulus  * 

0.40-0.50-0.60 

32-50-69 

Altocumulus 

0.37-0.50-0.69 

33-52-90 

Mimbostratus 

0.59-0.66-0.77 

31-49-62 

Cirrus/cirrostratus 

0. 12-0.31-0.62 

14-36-62 

;  Cumulonimbus  developing  or  mature 

0.75-0.76-0.92 

62-71-65 

1  Cumulonimbus  decaying 

0.66-0.73-0.92 

30-49-61 

|  Sunglints 

0. 10-0.24-0.40 

-56- 

,  Water 

0.03-0.06-0.01 

0.2-1. 5-4. 5 

Land 

» 

I 

0.05-0.09-0.16 

2. 5-5. 6-6.0 

Motes 

i 

Between  10  and  40  studied  objects  i 

n  every  class 

Sun  elevation  angle,  E,  between  30 

and  55  degrees 

*  seven  to  eight  oktas  (eighths)  of  sky  cover 


The  separation  of  sixteen  classes  in  a  four-dimensional 


intensity  space  (maximum  ana  minimum  values)  for  NOAA-6  ima 


elevation  angles  of  between  41.2  ana  47—  (after  Lilias,  1&&4) 


Cloud/surface 


Cb 

132-172 

100-150 

101-255 

163-255 

Cb 

132-172 

100-150 

101-255 

160-162 

Ns 

106-131 

60-130 

101-255 

160-255 

Cs 

48-105 

30-90 

101-255 

160-255 

Ci  over  water 

24-47 

15-30 

101-255 

121-1&2 

Ci  over  land 

24-47 

31-50 

101-255 

121-182 

Cu  congestus 

106-172 

80-150 

101-255 

140-159 

As 

46-80 

30-80 

101-255 

121-159 

Ac 

81-105 

65-100 

101-255 

140-159 

Sc 

81-131 

65-120 

101-255 

100-139 

Mist  over  land 

27-40 

38-60 

101-255 

60-120 

Mist  over  water 

24-40 

22-37 

101-255 

60-120 

St/fog 

41-60 

30-80 

101-255 

bO- 120 

hater 

0-23 

0-20 

101-255 

60-120 

Land 

0-26 

21-60 

101-255 

60-120 

Sunglints 

10-100 

10-100 

0-100 

60-120 

examples  given  the  diagnostic  value  of  the  classification  in  real-time 
meteorological  analysis  and  forecasting  appears  to  be  high. 

2.2  Statistical  Cloud  Retrieval 

In  either  clear  sky,  single  cloud  type  or  multiclouo  type 
situations,  projection  of  the  multispectral  signature  of  each 
satellite  image  pixel  onto  a  corresponding  multi-dimensional  histogram 
will  give  rise  to  the  presence  of  areas  of  relatively  high  pixel 
density.  These  regions,  dispersed  throughout  the  radiance  intensity 
space  are  associated  with  relatively  homogeneous  emitting  and 
reflecting  surfaces,  clouo  types,  oceans  and  land.  The  ability  to 
describe  and  to  discriminate  effectively  between  these  surfaces  or 
classes,  essential  for  any  meaningful  clouo  classification  scheme,  has 
prompted  the  invocation  of  statistica 1 ly-baseo  retrieval  algorithms 
which,  as  opposed  to  the  threshold  methods,  are  designed  to  treat 
large  groups  of  pixels  at  a  time,  corresponding  to  each  surface  type. 
The  methodology  employed  in  separating  spatial  radiance  patterns  into 
cloudy  and  clear  classes  can  be  thought  of  as  being  equivalent  to  the 
use  of  a  constant  threshold  and  clear  sky  radiance.  For  complex  cases 
this  would  appear  to  be  an  unattractive  concept,  however  as  wi.  11  be 
seen,  use  of  time  compositing  can  be  beneficial  for  statistical 
algorithms  as  well. 

The  fundamental  philosophy  of  statistical  cloud  analysis  is  the 
partitioning  of  multi,  dimensional  frequency  histograms  into 
representative  classes  for  which  three  recognised  routes  exist.  These 
are  (il  gaussian  histogram  analysis  (Simmer  et_  a  1.,  19&2;  Phulpi.n  e  t 
a  1. ,  1  963  );  (ii)  dynamic  clustering  (Lesbois  e_t  aj^,  19  621;  (i.ii.l 
spatial  coherence  method  (Coakley  and  Bretherton,  19b 2;  Coakley.  1963; 
Coakley  and  Baldwin,  19641.  As  the  different  method  descriptions  use 


the  terms 


"class"  and  "cluster"  synonymously  and  Interchangeably, 


this  convention  will  also  be  adoptea  here. 

2.2*1  Gaussian  Histogram  Analysis 

Simmer  e_t  a  1. ,  (1962)  along  with  Phulpin  e_t  aj^,  (  1  9 3 )  present 
examples  of  histogram  partition  by  way  of  fitting  Gaussian  (normal'' 
distribution  functions  to  one  ana  two-dimensional  frequency  histograms 
in  oraer  to  isolate  aistinct  clusters.  At  the  time  of  publication  the 
latter  method  haa  only  been  developed  ana  tested  for  cloua  layers 
overlying  oceanic  regions,  although  the  intention  was  to  aevelop  a 
similar  technique  for  use  over  land  surfaces.  The  asymmetric  Gaussian 
analysis  of  Simmer  et  a  1. , (  1 962)  does  not  suffer  from  this  limitation. 
Neither  method  at  present  takes  account  of  corrections  for  viewing 
geometry  or  atmospheric  effects.  Implicitly  built  into  these  analyses 
is  the  radiative  model  in  which,  as  previously  stated,  different 
surfaces  can  be  represented  as  peaks  in  the  histograms  by  virtue  of 
their  homogeneity  ana  consist  of  clusters  of  distinct  radiance  pairs. 

In  Phulpin  e_t  a  1«,  (1963)  the  AVhRR  data  array  was  first  aiviaea 
up  according  to  the  operational  grid  for  the  Tiros  operational 
vertical  sounder  (TGVS).  This  was  because  the  work  was  carried  out  in 
the  context  of  improving  vertical  temperature  profiles  in  clouay 
areas.  The  data  are  subsequently  calibrated  ana  the  visible  one¬ 
dimensional  and  two-dimensional  histograms  constructed.  Gaussian 
functions  with  fi.xea  standara  deviations  are  fitted  to  the  histogram 
peaks  to  identify  aistinct  sea  ana  cloua  clusters  which  must  contain 
at  least  a  threshold  minimum  number  of  pixels.  Closely  adjacent  peaks 
are  taken  to  represent  the  same  layer  and  are  merged  into  one  class, 
follow  i.ng  a  pre-aef  i.nea  set  of  rules.  Cloua  cover  estimation  is  by 


summing  up  all  the  pixels  in  the  clouay  clusters  ana  dividing  by  the 
total  number  of  array  pixels. 

Validation 

Validation  would  consist  of  comparison  with  visually  interpretea 
enhanced  images  by  experienced  cloua  analysts,  identifying  the  clouas 
by  their  texture  ana  spectral  features.  They  are  reckonea  to  be  able 
to  estimate  total  cloud  amount  to  within  an  accuracy  of  10%. 
Combinations  of  aata  from  channels  1  (0.5b  -  0.6b  yjm)  ana  4  (10.5  - 
11.5  |>m),  along  with  2  (0.73  -  1.10  y  m )  and  4  were  usea  in 
construction  of  the  histograms  with  the  result  that  the  near  infrarea 
channel  2  proved  significantly  better  than  channel  1  for  cloua 
analysis,  mainly  aue  to  the  lower  reflectance  of  the  sea,  thus 
allowing  for  better  spectral  resolution  of  cloua  free  regions  (see 
Figure  b).  Altogether  almost  90%  of  the  cases  were  in  close  agreement 
with  the  visual  interpretations,  with  the  most  persistent 
difficulties  arising  in  the  discrimination  between  small,  low-level 
cumulus  and  the  sea  surface,  a  problem  akin  to  that  causea  by  boundary 
layer  cumulus  in  the  threshold  method.  The  potential  role  of  channel 
3  (3.55  -  3.93  y-im),  although  consiaerea,  is  as  yet  unresolved. 

Ihe  idea  developed  by  Simmer  e_t  a  1. ,  (  1 9  b  2 )  is  that  peaks 

appearing  in  one-ai.mensiona  1  histograms  (both  visible  and  infrared'' 
can  be  represented  by  Gaussian  functions  which  can,  in  turn,  define 
clusters  in  two-dimensional  radiance  space.  The  standard  deviations 
of  the  Gaussian  curves  are  determined  from  the  plot  of  radiance  versus 
the  logarithm  of  infrared  frequency  and  the  first  cluster's  pixels 
redefined  on  the  basis  of  distance  from  the  peak  radiance.  This 
process  is  carried  out  first  for  the  Infrarea  histogram  ana 
subsequently,  whilst  only  using  the  pixels  assigned  to  the  first 
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Figure  8(a)  One  dimensional  frequency  histogram  of  A  N  =  N  -  N  * 
constructed  from  AVHRR  data  (channels  1,  0.,58  -  0.6fa  jim 

and  4,  10.5  -  11.5  ^4V1))  for  Europe  and  the  Atlantic  Ocean 

where  N  represents  the  cl<jud  cover  computed  from  the 
retrieval  algorithm  and  N  is  the  value  derived  from 
visual  nephanalysis  (after  Phulpin  et  al.,  1983). 
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cluster,  on  the  visible  histogram  so  as  to  define  the  first  two¬ 


dimensional  cluster.  The  whole  process  is  repeated  until  either  all 
image  pixels  are  assigned  to  a  cluster  or  a  set  number  of  clusters 
have  been  defined,  any  remaining  pixels  being  assigned  to  the  nearest 
(in  two-dimensional  space)  cluster.  Cloud  cover  fraction  is  estimated 
as  in  Phulpin  e_t  a  1.  ( 1  9 b 3 )  whilst  surface  ano  cloud  properties  are 
calculated  from  the  mean  radiance  over  a  particular  cluster  for  pixels 
residing  within  a  radius  of  one  standard  deviation  of  the  peak 
frequency.  Stable  statistics  can  be  achieved  by  operating  on  a 
subregion  of  at  least  300  image  pixels. 

This  latter  method  was  included  amongst  those  tested  in  the  1SCCP 
pilot  study.  Its  degree  of  success  was  judged  to  be  dependent  upon 
the  decision  criteria  ano  the  resulting  effectiveness  of  the  cluster 
definition.  It  performed  to  around  the  same  level  as  the  bispectral 
threshold  but  encountered  problems  in  regions  of  high  spatial 
inhomogeneity  and  in  the  case  of  a  time-varying  surface. 

2.2.2  Dynamic  Clustering 

A  significant  aspect  of  any  clouo  retrieval  is  that  it  is  carried 
out  in  the  most  objective  manner  possible.  The  examples  so  far  have 
involved  tasks  such  as  threshold  tuning  and  the  use  of  visual  image 
interpretation  which  invoke  a  degree  of  subjectivity.  Correspondingly 
the  most  objective  cluster  definition  should  give  rise  to  the  most 
objective  statistical  method.  The  dynamic  cluster  algorithm  (Desbois 
e  t  a  1. ,  19&2;  Desbois  and  Stze,  19b4a^  has  an  essentially  objective 
framework  and  has  been  modified  several  times  over  various  cloud 
distributions.  The  step  by  step  procedure  of  cluster  selection  can  be 
found  in  Desbois  et  a  l.,  (19&2''.  Here,  the  method  is  discussed  in  the 
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overall  context  of  classification  (Desbois  ana  Stee,  19&4a,b  ana  Sfeze 
et  al«,  1  9 b 4  ) . 

The  statistical  requirements  of  the  technique  involve  the  use  of 
several  thousanh  pixels  (termed  the  learning  set)  in  a  test  run. 
These  may  or  may  not  represent  the  image  pixels  on  which  the 
classification  will  be  performed.  One,  two  or  three-dimensional 
histograms  are  then  constructea  from  the  learning  set.  As  opposea  to 
Gaussian  techniques,  those  that  have  clusters  classified  from  hETEOSAT 
imagery  (as  in  Lesbois  e_t  aj^.,  19b2)  use  the  available  water  vapour 
channel  (5.7  -  7.1  jjm)  in  addition  to  the  visible  (0.4  -  1.1  yui)  and 
infrared  (10.5  -  12.5  jim),  resulting  in  the  use  of  three  insteaa  of 
two-dimensional  histograms.  Correct  alignment  of  the  image  channels  is 
neeaea  because  the  position  of  every  pixel  in  the  various  histograms 
will  determine  its  class.  Application  of  the  dynamic  cluster  method  to 
the  histograms  now  follows,  either  from  randomly  selected  pixel  groups 
(kernels)  or  from  pre-des ignatea  kernels.  Although  at  times  useful, 
the  pre-designated  kernels  are  to  be  avoidea  if  the  classification  is 
to  remain  wholly  objective.  The  end  product  of  the  procedure  is 
firstly,  the  number  of  classes  separated,  secondly  their  centres  of 
gravity  and  variances,  and  finally  the  percentage  of  pixels  from  the 
learning  set  in  each  class.  Classification  of  the  image  under  study  is 
achieved  by  assigning  each  pixel  to  the  class  it  lies  closest  to,  the 
distance  to  a  class  referring  to  the  sum  of  the  variance  of  the  class 
and  the  Euclidian  separation  between  the  image  pixel  and  the  centre  of 
gravity  of  a  class. 

Visual  examination  of  the  histograms  allows  a  qualitative 
assessment  of  the  cloud  systems  present  to  be  made  and  in  simple  cases 
the  classes  can  be  identified  directly  from  the  histograms.  Examples 
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of  one  and  two-oimensiona  1  histograms,  similar  to  those  given  in 
Desbois  et_  a  1.  (1^62''  have  been  derived  from  NGAA-7  AlhRR  data 
depicted  in  Figure  1.  The  one-dimensional  histograms  (Figures  9-12'' 
illustrate  peaks  of  varying  intensity  and  definition,  the  visible 
(channel  1)  data  tending  to  be  more  smoothed  than  the  infrared  (channel 
5),  resulting  in  easier  class  recognition  i.  n  the  latter.  Ihe 
separation  of  the  three  major  classes  present  in  the  sub-scene  termed 
here  frame  A  [  52°  -  57°]  (Figures  9  and  10),  (land  surface,  sea 
surface  and  low  stratocumulus)  is  more  easily  performed  on  the  infrared 
histogram.  Likewise  the  sub-scene  termed  frame  C  [  42°  -  ^47°]  (Figures 
11  ano  12),  the  very  broad  peak  lying  adjacent  to  the  land  surface  peak 
contains  both  middle  and  high  level  convective  clouds  which  can  only  be 
distinguished  on  the  infrared  histogram.  The  corresponding  (Figures  13 
ano  15)  two-dimensional  histograms  are  shown  along  with  the 

corresponding  results  for  the  sub-scene  termed  frame  b  147°  -  5 2°]  in 

* , 

Figure  1 U.  Once  again  several  classes  are  immediately  identifiable  by 
their  positions  in  the  radiance  intensity  space  whilst  those  less 
obvious  would  subsequently  be  deduced  from  dynamic  clustering.  broader 
peaks  in  the  one-dimensional  histograms  express  themselves  as  areas 
contained  within  but  which  are  themselves  devoid  of  contours  in  two- 
dimensions.  Projection  of  the  two-dimensional  plots  onto  a  three- 
dimensional  surface  (Figure  16)  provides  an  additional  aid  to  visual 
class  separation. 

The  major  drawbacks  concern  occasions  when  the  histograms  fail  to 
show  definite  peaks  but  instead  tall  off  gradually  from  a  single,  ill- 
defined  peak  (Uielicki  ano  Welch,  196b).  In  such  a  case  the  algorithm, 
whilst  still  separating  classes,  fails  to  draw  sharp  boundaries  between 
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Figure  S  One  dimensional  frequency  histogram  for  the  AVHRR 

visible  channel  1,  (0.5&  -  0.6b  ^m)  for  sub-scene,  frame 

A  (52°  -  57°)  of  data  shown  in  Figure  1. 
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Figure  1 1 


0 n e - d i m e n s i o n a  1  frequency  histogram  for  the  visible 
AVHRR  channel  1  (0.58  -  0.68  pm)  for  sub-scene,  frame  C 
(42°  -47°)  of  data  shown  in  Figure  1. 
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Figure  1  2 


One-dimensional  frequency  histogram  for  the  infrared 
AVHRK  channel  5  (11.5  -  12.5  >>m)  for  sub-scene,  frame  C 
(42°  -  47°)  of  data  shown  in  Figure  1. 
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re  22  Infrared/ visible  {AVHRR  channels  1,  0.56  -  0.68  um  and 

5,  11.5  -12.5  um)  bidimensiona 1  histogram  for  sub-scene 

frame  A  (5 2°  -  57°)  of  data  shown  in  Figure  1. 
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Figure 


I  n  f  ra  re  d  /  v  i  s  i  b  1  e  (channels  1,  0.56  -  0.66  yw  and  5, 

11.5  -  12.5  um)  b i d i me n s i o na 1  histogram  for  sub-scene, 

frame  B  (47°  -  52°)  of  data  shown  in  Figure  1. 
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Figure  1  5  Infrared/visible  (channels  1,  0.56  -  0.68  pm  and  5, 

11.5  -  12.5  pm)  b i d i m e n s  i  o n a  1  histogram  for  sub-scene, 
frame  C  (42°  -  47°)  of  data  shown  in  Figure  1. 


neighbouring  classes.  This  effect  was  nctea  in  the  ISCCF  pilot  study 
(R os sow  e t  a l « ,  1985)  ana  would  produce  analogous  problems  in  the 

Gaussian  analyses  (Section  2.2.1).  Secondly  the  dynamic  cluster  method 
requires  initial  choices  of  the  number  of  classes  and  the  number  of 
points  in  each  kernel.  These  values  can  have  a  profound  effect  upon 
the  number  of  final  classes  that  result  which  is  generally  too  large. 
Application 

Application  of  this  type  of  cloud  classification  to  individual 
image  segments  carries  with  it  the  inherent  risk  that  adjacent  segments 
might  not  produce  equivalent  classes,  depending  on  the  cloud  amount  ano 
distribution.  The  same  can  be  said  of  the  same  image  segment  examined 
at  various  time  intervals.  Should  this  occur  it  would  make  larger 
scale  classification  extremely  difficult.  Gne  solution  to  this  problem 
is  to  begin  with  learning  sets  of  a  sufficient  size  so  as  to  exceed  the 
minimum  number  of  points  required  for  a  representative  cloud  population 
ano  statistical  stability.  In  this  way  the  classes  found  at  larger 
scales  can  be  reproduced  at  smaller  scales.  Alternatively  time 
sampling  can  be  used  whereby  the  learning  set  is  derived  from  time 
series  of  an  image  segment,  each  image  taken  at  the  same  time  each  day 
for  several  successive  days.  The  resulting  cumulative  histograms  are 
then  likely  to  contain  distributions  of  most  types  of  clouo  and 
surface.  An  interesting  feature  of  this  sampling  would  be  the  large 
number  of  partially  cloudy  pixels  spread  out  through  the  histograms,  as 
would  be  those  pixels  depicting  semitransparent  cloud. 

The  original  method  was  tested  on  h&TEGSAT  iuiages  over  the 
tropical  Atlantic,  central  ano  West  Africa  (Lesbois  e_t  a  1. ,  1982  ). 

Although  no  corrections  for  atmospheric  effects  or  viewing  geometry 


were  made,  the  close  proximity  of  the  stuoy  region  to  the  equator 
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reduced  the  likelihood  of  the  latter  type  of  error.  Individual  image 
segments  used  for  the  initial  testing  for  statistical  stability 
included  bOGO  pixels  and  the  subsequent  full  study  involved  40,GGG 
pixels.  The  analysis  centered  on  an  area  of  convective  disturbances 
providing  a  varied  class  distribution.  The  retrieval  produced  up  to  o 
different  classes  (see  Figures  17  ano  16  and  Table  iJ)  which  were 
identified  by  referring  to  the  three  two-dimensional  histograms  ano  the 
rectified  images.  The  classes  were  found  to  correspond  well  with  the 
clouo  types  present,  encouragingly  discriminating  the  presence  of 
cirrus. 

In  addition  to  classification  Desbois  et  al»,  (19b£)  demonstrate 
a  method  of  determining  semitransparent  cloud  top  temperature  from  the 
infrared/water  vapour  histogram.  Results  were  encouraging  and  could 
possibly  be  used  as  correlative  data  in  future  applications  of  other 
retrieval  algorithms. 

Methods  for  Improving  Classification 

Since  the  original  testing  of  the  dynamic  clustering  technique 
reported  in  1 9 2: ,  much  effort  has  been  devoted  to  finding  means  of 
achieving  greater  representation  in  the  clouo  classes.  Some  new  ideas 
have  since  emerged:  (i)  use  of  space  and  time  sampling  to  obtain  a 
fuller  representation  of  the  clouo  distribution  (discussed  in  the 
previous  section)  (Desbois  and  Skze,  19&^a)  and  (ii)  comparison  of 
local  spatial  variances  to  separate  the  more  homogeneous  clouo 
clusters  (Size  and  Desbois  1  9  &  6 )  •  One  of  the  objectives  of  Desbois 
ano  Sfeze(  lS&4a,b),  was  to  perform  a  modified  retrieval  using  the 
recommendations  above  ano  to  compare  the  results  with  those  derived 
from  thresholds,  Ihe  first  modification  is  that  of  deriving  the 
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Classes  obtained  by  the  dynamic  cluster  algorithm  for 
the  same  analysis  region  as  in  Figure  17,  represented  on 
two-dimensional  histograms  IR-VIS  and  IR-WV  (after 
D  e  8  b  o  i  s  e_t  ji_lj_,  1  98  2  ) . 
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1  7  and  1  8.  Percentage  coverage  o f  each  class  along  with  the  centre  o f 
gravity  and  variance  of  each  class  are  Ktven.  Ihe  units  are  numerical 
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surface  radiances  from  time  series  of  images*  Each  image  is  extracted 
at  the  same  time  each  day  in  order  to  avoid  the  requirements  for 
corrections  for  variations  in  solar  zenith  angle.  Also  the  use  of  raw 
METEOSAT  radiometric  counts  for  analysis  avoids  the  use  of  calibration 
factors  and  angular  and  atmospheric  corrections,  all  of  which  vary  in 
accuracy.  These  stipulations  amount  to  the  requirement  that 
observation  conditions  remain  approximately  constant  throughout  the 
period  of  image  gathering.  Size  and  Desbois  (  1  9  b  b )  discuss 
compositing  techniques  which  are  made  rather  complicated  by  the  lack 
of  correlation  between  minimum  reflectance  and  maximum  brightness 
temperature  from  pixel  to  pixel,  a  feature  that  often  arises  in  the 
presence  of  very  thin  cirrus  cloud.  It  is  worth  noting  that 
information  about  daily  cloud  distributions  ana  temperature  maxima  can 
also  be  ascertained  from  the  inf rared/ v isib le  histograms  when  certain 
cloud  types  are  fairly  clearly  defined,  such  as  cirrus,  low  level  and 
multilayered  clouds. 

The  use  of  spatial  variances  is  designed  to  discriminate 
spectrally  similar  regions  as  well  as  provide  information  on  partially 
cloud  covered  pixels  (Coakley  and  Bretherton,  19bB^  (i.e.  those 

regions  where  classes  tend  to  overlap  in  two-dimensional  histograms). 
Clouds  (surfaces)  of  variable  optical  depth  but  uni. form  emissivity 
will  have  a  high  visible  variance  but  low  infrared  variance  such  as 
low  cumuliform  clouds.  The  reverse  applies  to  clouds  of  uniform 
optical  depth  ana  variable  emissivity  such  as  semi-transparent  cirrus. 
The  use  of  v isib le/ inf rared  variance  histograms  is  complementary  to 
the  radiance  histograms  which  suggests  that  a  four  (five''  parameter 
classification  lh/IRva„/VXS/VISvav,  (lh/Ihva„/V  IS/V  ISva„/V»  V  )  may  help 
to  separate  further  those  cloud  classes  already  broadly  defined  in  the 


original  classification.  Such  a  rr.oaifiea  dynamic  cluster  analysis 
using  four  parameters  is  described  in  Sfeze  ana  Lesbois  (15&6''.  In  this 
version  the  use  of  local  variances  improves  discrimination  of  cirrus 
ana  improves  the  separation  of  totally  clouay  classes  from  partially 
covered  pixels. 

Ihe  problem  of  treating  partially  clouay  pixels  i.  s  well 
established.  A  novel  approach  to  solving  this  problem  has  been 
proposea  by  Arking  ana  Chilas  ( 1 9  b  5 )  *  Their  methoa  is  aesignea  to 
treat  pixels  lying  outsiae  aefined  clusters  ana  determine  their  cloua 
cover  fraction  whilst  requiring  that  their  other  cloua  parameters, 
optical  depth  ( S)  and  cloud  top  temperature  (Tc),  remain  as  consistent 
as  possible  with  those  of  their  nearest  cluster. 

The  algorithm  requires  information  from  three  AVhhh  channels  (1. 
0.5&  to  0.6 b  yjm;  35  3«55  to  3*93  um  and  10.5  to  11.5  ^m)  ana  in  the 
process  retrieves  four  parameters,  f,S,Tc  and  a  microphysical  moael 
index  m  describing  the  distribution  of  size,  shape  and  refractive 
index  of  the  cloua  particles  present  in  each  pixel.  Six  possible 
values  of  m  have  been  computea  (Table  5).  The  basic  raaiative  model 
assumes  only  one  cloud  type  per  pixel  ana  uniform  horizontal  ana 
vertical  temperature  and  optical  properties  within  each  pixel.  Ihe 
retrieval  (£,tc  and  m  are  determined  once  f  is  known)  uses  equations 
relating  the  three  channel  raaiances  measured  for  each  pixel  to  the 
scene  parameters  within  that  pixel.  A  second  radiative  transfer  model 
is  thus  built  into  the  algorithm  which  i.  s  essentially  a  four-stage 
process,  begi.  nni.  ng  with  the  determination  of  the  clusters  from 
histograms  ana  subsequently  the  surface  parameters  for  each  channel. 
The  microphysical  moael  1  n a e x  appropriate  to  each  pi x e  1  i. s  founa  by 
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Table  5.  Nicrophysical  models  (NiPhs)  and  their  optical 


parameters 


MPM 

index 

(after 

Arkinft 

and  Chiles, 

1965). 

Mode 

radius 

(ym) 

TIftOS-N 
Particle  channel 
phase  (^im) 

Extinction 

coefficient 

(cnr/particle) 

Particle 

albedo 

Asymmetry 

factor 

0*73 

1.66E-6 

1.000 

0.645 

1 

4 

liquid 

3.7 

2.01E-6 

0.937 

0.753 

11 

0.66E-6 

0.5&3 

0.656 

0.73 

1.68E-6 

1.000 

0.663 

2 

4 

ice 

3.7 

1.96E-6 

0.884 

0.751 

11 

1.16E-6 

0.305 

0.643 

0.73 

b.55E-6 

1.000 

0.663 

3 

8 

liquid 

3.7 

7.18E-6 

0.676 

0.620 

11 

5.93E-6 

0.494 

0.536 

0.73 

6.55E-6 

1.000 

0.671 

4 

8 

ice 

3.7 

7.71E-6 

0.795 

0.637 

11 

5.97E-6 

0.411 

0.527 

0.73 

2.57E-5 

1.000 

0.873 

5 

16 

liquid 

3.7 

2.73E-5 

0.601 

0.873 

11 

2.83E-5 

0.514 

0.566 

0.73 

1.02E-3 

1.000 

0.887 

3.7 

1.06E-3 

0.606 

0.936 

comparing  the  observed  channel  3  radiance  to  that  calculated  for  each 
'  ’■'dex  using  model-derived  values  of  £  and  Tc  (Figure  15). 

Values  of  optical  depth  and  clouo  top  temperature  can  be 
determined  for  all  ’cluster-pixels’  by  application  of  the  radiative 
transfer  equations  that  effectively  map  visible  reflectance  and 
infrared  brightness  temperature  into  £  ano  1 c  assuming  f  and  m  are 
known  (f=l  for  defined  clusters).  Recovery  of  whole  image  clouo 
statistics  cannot  be  achieved  however,  until  the  values  of  f  for 
partially  clouay  pixels  (i.e.  those  lying  outside  the  clusters)  have 
been  determined.  Maximal  clustering  is  now  invoked  as  follows. 
Figure  20  illustrates  curves  constructed  from  radiative  transfer 
models  describing  the  variation  of  £  and  Tc  with  clouo  cover  fraction 
for  a  given  visible  reflectance,  infrared  brightness  temperature  and 
microphysical  model  index.  The  position  of  each  pixel  is  allowed  to 
vary  along  its  respective  curve  until  it  lies  closest  to  one  of  the 
clusters,  in  this  way  maximising  the  pixel  concentration  around  the 
clusters.  The  value  of  f  that  corresponds  to  the  new  position  of  the 
pixel  along  its  curve  is  set  as  the  retrieved  value  of  f,  permitting  <£ 
ano  lc  to  be  determined  (Figure  21). 

The  algorithm  was  tested  on  one  of  the  1SCCP  pilot  study  test 
regions,  large  enough  to  ensure  stable  statistics  and  a  satisfactory 
level  of  consistency  reached  in  a  comparison  with  manually  interpreted 
images.  Use  of  the  improved  clustering  method  described  would  be  a 
I  logical  improvement  ano  a  direct  comparison  to  results  from  the  other 

clustering  methods  could  serve  to  widen  knowledge  of  this  technique. 
2.3  Spatial  Coherence  Algorithms 

The  major  impetus  for  the  study  of  local  spatial  variances  (ana 
their  subsequent  inclusion  into  dynamic  cluster  algorithms!  was 
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Figure  20 


Possible  values  of  the  cloud  parameters  S  and  Tc 
associated  with  a  given  11  micrometers  brightness 
temperature  (T^=260K)  and  two  visible  reflectances 
(Rv=0.30,  0.55)  (after  Arking  and  Childs,  1985). 


provided  by  Coakley  ana  bretherton  (15b 2)  who  aiscoverea  that  local 
variances,  as  applied  to  simple  layered  cloua  systems  coula  retrieve 
the  raaiances  corresponding  to  the  cloua  free  and  c loud-coverea  fielas 
of  view,  as  well  as  those  emanating  from  partially  fillea  fielas  of 
view.  The  fundamental  assumptions  underlying  their  spatial  coherence 
procedure  are  that  a)  the  clouds  are  situatea  in  aistinct  layers  with 
each  layer  possessing  a  temperature  appropriate  to  its  altitude,  ana 
b)  all  the  clouds  in  a  layer  emit  radiation  that  is  characteristic  of 
the  layer,  the  layers  being  optically  thick  at  the  observea  wavelength 
with  zero  transmissivities.  This  being  the  case,  when  several  cloua 
layers  occur,  each  layer  having  at  least  several  pixels  wiath  of  non- 
overlappea  cloua,  the  methoa  permits  the  estimation  of  the  radiances 
from  the  completely  clouay  regions  of  each  layer.  by  allowing 
detection  only  in  the  infrarea  channel  the  retrieval  can  be  usea  to 
proauce  aiurnal  analyses. 

To  illustrate  the  simplest  case  the  local  standara  aeviation  of 
11  micrometer  brightness  temperatures  is  plotted  in  Figure  22  as  a 
function  of  local  mean  radiating  temperature  for  6x6  pixel  arrays 
(about  32km  x  32km)  of  AVHRR  global  area  coverage  oata  for  an  oceanic 
region  underlying  a  simple,  homogeneous  cloua  layer.  The  result  is  an 
arch-like  structure,  typical  of  such  plots,  comprising  two  distinct 
clusters  of  low  stanaara  aeviation  points  separated  by  a  dispersion  of 
higher  standard  aeviation  points.  The  ’warm'  cluster  centered  on  293k 
+ 1  k  is  attributed  to  homogeneous  clear  sky  fielas  of  view  (ocean 
surface)  whilst  the  colder  cluster  centerea  around  2b 3.5k  +  1.5k 

corresponds  to  raaiances  emittea  from  completely  clouay  pixels  ana  the 
height  of  the  cloud  layer  can  be  founa  from  the  difference  in  radiating 
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Figure  22  Local  means  and  local  standard  deviations  for  the  3«7 

micrometers  and  11  micrometers  channels  for  a  simple  one¬ 
layered  cloud  system  over  the  Pacific  Ocean  (after  Molnar 
and  Coakley,  1 9 fcs 5  ) • 


temperature.  The  points  residing  in  the  boay  of  the  arch  are 
attributed  to  partially  filled  pixels  (with  f  decreasing  as  we  increase 
the  radiating  temperature  from  2b 3.5K)  as  these  pixels  should  not  be 
expected  to  possess  the  same  degree  of  local  coherence  as  those  from 
clear  or  cloudy  regions.  hence  since  the  presence  of  a  layer  can  only 
be  detected  when  adjacent  fields  of  view  are  completely  filled  this 
retrieval  method  is  successful  only  when  the  cloud  cover  is  uniform 
over  areas  considerably  larger  than  the  resolution  of  the  data.  If 
this  crieterion  is  not  met,  the  resulting  plot  could  not  be  used  to 
determine  the  cloud’s  radiating  temperature.  Those  cloud  types 
exhibiting  considerable  spatial  structure  such  as  developing  cumulus 
congestus  and  cirrus  are  therefore  less  likely  to  be  detected  as  a 
layer. 

Estimation  of  fractional  cloud  cover  ( AQ )  for  regions  comprising 
single  layer  systems  is  given  by 

Ac  =  (I  -  IS)/(IC  -  Is)  (1) 

where  I  is  the  mean  radiance,  defined  by 

1  =  <1  -  M  Is  +  Me  (2  > 

where  ls,  the  clear  sky  radiance,  and  Ic,  the  cloudy  radiance  are  both 
derived  from  the  position  of  the  ’’feet"  of  the  arches.  The  uncertainty 
in  Ac  can  be  estimated  once  the  uncertainties  in  Is  and  Ic  are  known 
since  they  are  proportional  to  the  breadth  of  the  arch  feet.  Any 
uncertainty  in  I  is  the  result  of  instrument  noise  and  is  of  negligible 
magnitude.  Note  that  the  method,  applied  to  clear  sky  conditions  over 
oceanic  regions,  will  provide  a  useful  guide  to  sea  surface  temperature 
(by  applying  the  inverse  Planck  function  to  the  clear  sky  radiance 
value),  although  with  greater  uncertainty  compared  to  values  obtained 
from  sounding. 
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Because  persistent  single-layered  cloud  regimes  (e.g.  marine 
stratocumu lus  near  coastal  South-West  Africa)  are  fairly  scarce  ana 
the  majority  of  cloud  systems  around  the  globe  comprise  two  or  more 
distinct  layers,  the  ideas  of  spatial  coherence  were  extendea  to 
multilayered  systems  (Coakley,  1 9  6  3 )«  When  two  layers  appear  in  a 
scene  the  points  in  the  bodies  of  the  resulting  arches  aenote 
contributions  of  cloud  cover  from  both  layers  as  in  Figure  23. 
Specifically  in  Figure  23(a)  (Figure  23(b))  the  regions  of  locally 
uniform  emission  are  observea  at  approximately  bb,  6b  ana  37  mW  v^~eL 
s  r“  ^  cm  (0*30,  0.14  and  0.03  mW  m“^  sr"^  cm).  The  arch  foot  at  bb 
(0.30)  mW  m  sr"  cm  is  interpreted  as  being  due  to  radiation  from 
the  cloud  free  oceanic  background  whereas  the  two  remaining  feet  are 
interpretea  as  lower  and  upper  level  cloud  layers  respectively.  All 
the  pixels  not  yet  accounted  for  are  designated  as  partially  coverea 
with  an  as  yet  unknown  combination  of  clouds  from  the  two  layers.  If 
the  fractional  contribution  of  each  layer  is  known  then  the  total 
fractional  cover  can  be  found.  Figure  24  shows  a  spatial  coherence 
plot  for  the  sub-scene  t  rame  C  [42  °  -  47°]  of  the  NCAA  AVhRfi  data  from 
Figure  1,  aepicting  some  uncertainty  in  the  aefinition  of  the  arches. 
Three  arch  feet  are  discernible,  that  centerea  on  3b0  counts  referring 
to  the  land  surface  background  whilst  the  major  cloud  decks  present 
give  rise  to  the  feet  centered  at  570  and  7b0  counts.  The  breaath  of 
the  arch  feet  and  the  large  number  of  pixels  within  the  (il 1-def ineal 
arches  suggest  that  in  this  case  a  poor  result  for  Ac  woula  be 
obtained.  This  illustrates  some  of  the  problems  of  performing  spatial 
coherence  retrievals  against  (i)  a  non-uniform  land  surface  ana  (ii)  a 
rather  complex  cloud  pattern,  Coakley  (1963)  usea  simultaneous 
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observations  in  the  3«7  micrometers  window  (channel  3  AVhhfO  ana 
the  11  micrometers  window  (channel  A  AVHRft),  (these  channels  being 
linearly  independent  from  each  other)  to  obtain  the  non-overlapping 
fractional  cover  attributed  to  each  of  the  two  layers  by  solving  pairs 
of  simultaneous  linear  equations  such  as  those  below.  for  observations 
in  channels  I  and  J,  and  for  two  cloud  layers  1  (the  lowest  layerl  ana 
2  (the  upper  layer) 


where  l1  is  the  radiance  received  in  channel  i;  Ac1  is  the  fractional 
area  of  a  pixel  covered  by  clouds  from  layer  1  that  are  not  overlapped 


by  layer  2  clouds;  Ac2  is  the  fractional  pixel  area  covered  by  only 
layer  2  clouds;  Is1  is  the  clear-sky  radiance  in  channel  i  and£-1iIc-)1 
represents  the  radiance  observed  for  pixels  cloud-coverea  by  layer  1 
with  no  overlapping  layer  2  clouds,  1  being  the  emissivity  of  the 
optically  thick  layer  1  cloud  in  channel  i.  Similarly  G.2  llc2  1  is  the 
radiance  from  pixels  completely  covered  by  layer  2  clouds.  Whilst 
and  Ic1,J  values  are  not  known  separately  the  values  of  their  product 


terms  in  Equations  (3)  and  (4)  are  derived  from  the  feet  of  the  arches, 
as  in  Figures  23  and  24.  Solution  of  Equations  (3)  ana  (4)  for  Acl  ana 
AC2  is  straightforward  in  an  ideal  situation. 

However,  with  the  likely  uncertainties  in  I  I  ^  and  lc2  noted 
Coakley  ( 1 9 6 3 )  has  aetermined  a  stability  conaition  for  the  equations 


which  indicates  when  they  cease  to  be  li.  nearly  independent. 


Satisfaction  of  this  conaition  entails  that  whatever  applies  to  two 
layers  can  apply  equally  to  three.  However,  the  number  of  linearly 
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independent  and  widely  spaced  infrared  channels  is  a  limiting  factor, 
indeed  the  use  of  the  3  •  7  micrometers  and  11  micrometers  channels 
together  permits  only  nighttime  observations  oue  to  daytime  reflection 
of  solar  radiation  by  low-level  water  clouds  in  channel  3»  The  window 
at  2.3  micrometers  is  as  yet  unexplored  but  may  be  studied  in  future 
from  Landsat  measurements.  Of  high  importance  is  the  need  to  test  the 
validity  of  the  underlying  assumptions  prior  to  testing.  Coakley's 
major  conclusion  was  that  cloud  systems  forming  arches  from  a  spatial 
coherence  analysis  almost  always  satisfy  both  layer  and  optical 
thickness  assumptions.  Data  presented  in  Coakley  (19&3)  were  extracted 
from  tropical  and  mid-latitude  Pacific  regions  and  excluded  continental 
areas  where  clear  sky  properties  are  more  variable.  The  results  showed 
promise  but  would  benefit  from  comparison  with  in-situ  data. 

The  suitability  of  the  spatial  coherence  technique  to  maritime 
environments  was  further  underlined  by  Coakley  ano  Baldwin  (1964)  who 
adopted  a  strategy  of  analysing  smaller  (60  km  x  60  km)  subregions  in 
preference  to  250  km  x  250  km  areas,  compositing  the  subregion 
statistics  to  provide  larger  scale  results.  Apart  from  reducing  the 
data  input  volume  by  more  than  an  order  of  magnitude  the  cloud 
patterns  encountered  in  the  subregions  were  almost  always  one  or  two- 
layered  ano  were  thus  conducive  to  analysis  by  spatial  coherence. 
Complex  systems  such  as  frontal  clouds  and  large  convective  systems 
(those  for  which  spatial  coherence  cannot  provide  real  i  Stic 
retrievals)  were  treated  with  a  default  estimate  [see  Coakley  ano 
Baldwin  (1964)  for  full  description].  This,  by  definition, 
incorporated  greater  uncertainty  than  the  spati.  a  1  coherence 
retrievals.  For  the  area  studied,  those  oceanic  regions  which  were 
frequented  by  the  ’awkward’  systems  suffered  from  larger  uncertainties 


in  mean  cloudiness  as  comparea  with  those  areas  dominated  by  layered 
cloud  formations. 

Validation 

This  approach  suggests  itself  as  a  possible  ISCCP  maritime 
retrieval  algorithm  (at  least  in  areas  where  the  climatological  cloud 
patterns  indicate  its  potential  applicability),  but  the  results  for 
landmasses  (although  yet  to  be  analysed)  are  unlikely  to  be  so 
promising,  particularly  in  mid- latitudes.  An  alternative  channel  to 
3*7  micrometers  needs  to  be  found  for  daytime  retrievals  lone  recently 
proposed  alternative  is  the  snow  /  cloud  discrimination  channel  centred 
at  1.6  micrometers]  and  corrections  for  viewing  geometry  have  not  been 
considered . 

The  most  recent  development  in  this  area  is  the  work  of  foolnar 
and  Coakley  (19&5)  who  assembled  data  sets  for  single  and  two-layereo 
cloud  systems  for  which  the  spatial  coherence  technique  provided 
reliable  results  ana  which  were  statistically  independent  of  one 
another.  bsing  these  examples  they  were  able  to  demonstrate 
similarities  in  cloud  cover  patterns  for  large  and  small  sub-scenes. 
This  discovery  may  be  of  interest  to  climate  modellers  who  wish  to 
represent  the  effect  of  sub-grid  scale  clouds.  The  analysis  was 
subsequently  used  to  design  a  statistically-based  algorithm  which, 
using  simple  mathematical  models  is  intended  to  fit  the  correct  value 
of  mean  cloudiness  over  all  the  partially  cloudy  pixels.  Initial 
results  were  comparea  to  those  derived  from  spatial  coherence  with  good 
consistency  achieved.  Total  cloud  amount  was  slightly  underestimated 
in  double-layered  cases  and  similar  errors  were  predicted  for  cloud 
systems  of  greater  complexity.  because  the  method  had  no  fundamental 


dependence  on  spatial  coherence  assumptions  holnar  ana  Coakley  ( 1 9 fa 5 ) 
emphasisea  its  potential  value  in  "difficult"  areas*  Results  of 
further  tests  are  now  required*  Overall  the  status  of  spatial 
coherence  methods  for  cloua  retrieval  is  one  of  limited  success  so  far 
but  clearly  worthy  for  future  development* 

2.4  Radiative  Transfer  -  type  Algorithms 

This,  the  third  major  class  of  retrieval  algorithms  employs  an 
entirely  different  methodology  from  the  retrieval  techniques  aescribed 
so  far.  Here  the  cloud  parameters  are  determinea  from  the  fitting  of 
atmospheric  radiative  transfer  models  to  the  observea  radiances.  Ihe 
0ioaels  aiffer  in  their  complexity  according  to  the  number  of  wavelength 
channels  usea  and  most  employ  some  auxiliary  data  to  improve  the 
accuracy  of  retrieval  (Reynolas  and  Vender  Haar,  1977;  Platt,  1 9 fa 3 ; 
Rossow  et  a  1.,  19b5)* 

hethods  of  calibrating  absolutely  infrared  radiances  from  several 
spectral  channels  against  appropriate  models  in  order  to  infer  cloud 
properties  has  been  put  forwara  by  Chahine  (  1962)  as  a  co-product  of 
vertical  temperature  profile  and  humidity  distribution  retrievals.  In 
this  method  the  technique  is  to  obtain  first  the  cloua  free  radiances 
by  comparing  the  radiances  measured  in  adjacent  areas,  assuming 
variations  in  the  cloud  properties  within  the  comparison  region.  The 
model  employea  describes  both  cloud  radiative  properties  ana  cloud 
vertical  distribution  ana  is  usea  in  the  comparison  of  the  observea 
radiances  to  the  estimated  cloud  free  values  to  extract  the  cloud 


properties. 

As  the  number  of  spectral  channels  used  increases  so  increasingly 
complex  cloud  models  can  be  employed  and  it  should  be  possible,  with 
sufficient  channels,  to  model  effectively  multilayered  cloud  systems 
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and  increase  the  accuracy  of  retrieval  in  such  cases*  At  present, 
however,  the  procedure  utilises  only  three  channels  ana  can  aerive  f 
ana  cloud  top  height  for  simple  (i.e.  single  cloud  layer)  models. 
Tests  so  far  inaicate  gooa  agreement  in  the  case  of  optically  thick 
cloua  with  large  deviations  encounterea  for  thin  high  clouds.  This 
difference  in  results  for  low  ana  high  clouds  has  hinderea  the  progress 
of  this  type  of  radiative  transfer  algorithm.  Future  improvements  rest 
on  the  inclusion  of  extra  spectral  bands  ana  model  refinements. 

The  bispectral  technique  aeveloped  by  Reynolas  ana  Vonder  haar 
(1977)  examines  each  pixel’s  net  radiation  balance,  relating  the 
radiances  received  in  both  visible  and  infrared  channels  to  the  sum  of 
the  contributions  from  clear  and  cloudy  regions  within  the  field  of 
view.  Input  parameters  to  the  model  include  surface  temperature, 
visible  albedos,  infrared  emissivities  and  visible  bidirectional 
reflectances  which  are  usea  to  calculate  the  fractional  radiance 
contributions  for  an  idealised  single  layer  cloud.  The  secona  step 
uses  the  radiances  to  calculate  the  cloua  fraction  ana  cloud  top 
temperature  for  each  pixel.  Vertical  temperature  profile  retrievals 
can  be  used  to  equate  Tc  with  cloud  top  height.  This  method  allows 
correction  for  viewing  geometry  and  varying  surface  properties, 
important  features  which  the  previously  aescribed  algorithms  ao  not 
incorporate.  The  test  site  used  haa  the  advantage  of  a  high  density  of 
surface  observers  ana  in  a  comparison  stuay  very  good  agreement  was 
reached  with  them  with  the  exception  of  cirrus  where  local  variations 
in  optical  depth  were  misinterpreted  as  variations  in  clouatop  height 
ana  f.  The  model  had  assoc iatea  a  constant  emissivity  for  cirrus 
clouds  but  subsequently  a  new  algorithm  was  devised  for  cirrus  clouas 


in  which  an  attempt  was  made  to  calculate  emissivity  as  a  direct 
function  of  visible  albedo#  This  greatly  inipro veo  cirrus  retrievals# 
The  other  major  ’problem*  cloud  type,  that  of  bounaary  layer  cumulus 
was  well  retrieved  by  this  technique. 

A  ’dynamic'  radiation  budget  approach,  as  opposed  to  the  'static' 
radiation  balance  described  above,  has  been  proposed  by  Smith  ano 
Vonoer  haar  (19&3)«  The  improvement  here  concerns  the  determination  of 
the  cloud  ano  surface  radiative  properties  from  the  satellite  data 
itself  instead  of  from  'static*  model  comparisons.  Much  of  the 
methodology  appears  in  previously  discussed  algorithms#  Firstly,  the 
surface  radiative  properties  are  deduced  from  time-period  composites 
(15  days)  of  minimum  reflectance  and  maximum  brightness  temperature  for 
each  pixel  location#  The  composites  must  then  be  spatially  averaged  if 
mean  surface  radiative  properties  are  to  be  determined  for  small  sub- 
regions.  Application  of  an  appropriate  clustering  technique  (Everitt, 
I960)  to  bispectral  radiance  histograms  generates  clusters  representing 
the  various  cloud  types  present,  each  cluster  center  defining  a 
particular  visible  albedo  and  cloud  top  temperature.  Climatological 
temperature  profiles  are  used  to  equate  each  Tc  with  a  clouo  top 
height# 

A  similar  radiative  model  to  that  used  in  Reynolds  and  Vonder 
Haar  (  1  977)  is  used  to  estimate  cloud  cover  fraction,  f.  It  is 
assumed  that  the  measured  mean  radiances  are  again  a  function  of  the 
combined  clear  and  cloudy  radiance  terms.  The  problems  of  variable 
infrared  emissivities  resulting  from  varying  optical  depth  are 
overcome  as  previously  by  the  inclusion  of  an  emissi vity-ref lectance 
function  whilst  corrective  terms  for  atmospheric  attenuation  in  the 
infrared  (mainly  cue  to  water  vapour  absorption  1  and  viewing  geometry 


are  also  included.  In  order  to  perform  on  a  diurnal  basis  the 
algorithm  employs  an  infrared  threshold  during  nighttime,  defined  so 
as  to  reduce  the  mean  square  difference  in  retrievea  daytime  cloud 
amount  between  itself  and  the  dynamic  bispectral  method,  thus  assuming 
that  the  most  suitable  threshold  can  be  used  irrespective  of  the  time 
of  day.  The  method  achieved  a  satisfactory  level  of  success  in 
preliminary  tests  with  the  exception  of  thin  cirrus.  Its  merit 

relative  to  its  predecessor  is  as  yet  unclear  and  a  direct  comparison 
would  be  a  useful  exercise. 

Occasionally  radiative  transfer  methods  include  threshold  or 
statistical  elements  in  their  framework  (Rossow  et  a  1»,  1  985).  A 

surface/atmosphere/single  cloud  layer  model  is  utilised  here  and  model 
radiances  calculated  as  a  function  of  the  surface  properties,  cloud 
properties,  vertical  profiles  of  temperature,  humidity  and  ozone 
abundance,  along  with  viewing  geometry. 

The  main  points  of  the  model  are  described  as  follows.  The 
atmosphere  is  composed  of  a  Rayleigh  scattering  gas  residing 
underneath  and  above  a  single  layer  of  cloud.  Its  vertical 
temperature  profile  is  determined  from  climatological  daily  average 
values.  The  cloud  layer  is  assumed  to  be  plane-parallel  comprising 
water  spheres  whose  effective  mean  radius  is  10  micrometers.  Visible 
reflectivity,  transmissivity  and  infrared  eraittance  are  all 
characteristic  of  such  a  surface.  Incoming  solar  radiation  is 
partially  attenuated  by  atmospheric  ozone  whose  distribution  is 
provided  from  a  seasonal  zonal  mean  climatology  (see  Hilsenrath  and 
Schlesinger,  1981).  The  land  surface  is  depicted  as  an  isotropic 
reflector  with  a  reflectivity  that  varies  with  position  whilst  all  the 
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surface  is  assigned  a  unit  emissivity  in  the  infrared.  The  ocean 
surface  likewise  has  unit  infrared  emissivity  whilst  its  reflectivity 
is  governed  by  Fresnel  reflection  coefficients  (see  Cox  and  hunk, 
1956).  Finally  weighted  reflectivities  are  introduced  to  any  region 
experiencing  snow  cover  (this  applied  to  one  of  the  test  regions  when 
the  model  was  employed  in  the  ISCCP  pilot  study). 

After  comparing  the  model  radiances  with  the  observed  values, 
simultaneous  infrared  and  visible  thresholds  (infrared  only  at  night) 
are  applied  to  the  deduced  cloud  properties,  the  criterion  being  that 
both  must  be  exceeded  in  order  to  count  a  pixel  as  cloudy.  Since  the 
thresholds  are  applied  to  cloud  properties  which  result  directly  from 
observed  versus  model  radiance  comparisons,  they  do  not  resemble  fixed 
radiance  increments  like  the  previous  thresholds.  The  results  from  the 
testing  of  this  algorithm  in  the  ISCCP  pilot  study  were  nominated  by 
its  repeated  failure  to  detect  low  lying  clouds  whilst  remaining  clouds 
were  retrieved  satisfactorily. 

A  third  bispectral  technique  involving  the  comparison  of 


VH 


ft 


3$ 


bispectral  curves  of  visible  albedo  versus  infrared  brightness 
temperature  from  "real"  clouds  to  model  curves  for  extraction  of  cloud 
fraction  and  optical  depth  is  given  in  Platt  (19&3)*  hodel  curves  have 
been  constructed  for  idealised  single  and  multilayered  clouns  and 
compared  to  bispectral  histograms  determined  from  observed  satellite 
data. 


Each  idealised  cloud  layer  is  characterised  by  an  optical  depth 
that,  although  constant  within  an  individual  pixel,  may  vary  between 
adjacent  pixels.  The  value  of  f  is  permittee  to  vary  from  one  pixel 
to  the  next  but  the  cloud  top  within  a  layer  is  set  at  a  constant 
height.  The  actual  cloud  is  defined  as  a  homogeneous  isothermal 
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region  whilst  scattering  effects  are  limited  to  incoming  solar 
radiation  which  the  layer  scatters  isotropically.  The  underlying 
surface  has  unit  infrared  emissivity,  a  visible  albedo  resembling  that 
of  an  ocean  surface  and  is  assumed  to  be  isotropic  as  well  as 
homogeneous.  Equations  for  the  visible  and  infrared  radiances  from 
'model'  clouas  that  might  be  received  by  a  satellite,  are  then 
formulated.  These  include  albedo,  emissivity,  pixel  cloud  fraction 
and  optical  depth  terms  ana  provide  the  theoretical  basis  upon  which 
the  curves  are  constructed  (Figure  25(a)).  The  parameters  used  in  the 
single  layer  calculations  are  shown  in  Table  6. 

The  curves  depict  albedo/brightness  temperature  relationships  for 
clouds  forming  (i)  a  continuous  layer  with  varying  optical  depth  and 
(ii)  a  broken  layer  with  uniform  optical  depth.  In  this  way  changes 
in  cloud  amount  and  optical  depth  can  be  related  to  changes  in  albeao 
ana  temperature.  Curves  1  and  2  could  be  said  to  depict  cruaely,  in 
the  high  cloud  case,  cirrus  and  cirrocumulus  clouds  respectively 
whilst  in  the  low  cloud  case  they  might  characterise  stratus  ana 
s tratocumu 1  us  respectively.  The  sensitivity  of  the  model  to  cloud 
detection  will  depend  on  the  albedo  and  temperature  resolution  of  the 
observations  and  on  the  data  array  having  a  smaller  dispersion  than 
the  distance  between  curves  1  and  2.  There  is,  of  course,  the 
additional  possibility  of  further  model  refinement.  Ihe  theory 
discussed  by  Platt  (1963)  also  permits  the  calculation  of  cloud 
fraction  and  emissivity  for  each  pixel  within  a  single  cloud  layer. 

The  method  has  been  extended  to  account  for  variations  of  albeao 
with  solar  zenith  angle,  variations  of  brightness  temperature  with 
satellite  viewing  angle  and  cases  where  multiple  cloud  layers  overlap 
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Figure  25(c)  Bispectral  curves  of  albedo, 


versus  brightness 


Figure  23(d) 


temperature,  Te,  for  two  overlapping  cloud  layers.  The 
layers  are  broken  but  cloudy  regions  have  uniform  optical 
depth  (after  Platt,  1  983). 

Bispectral  curves  of  albedo,  **<  ,  versus  brightness 
temperature,  Tg,  for  two  overlapping  cloud  layers  in 
which  the  top  layer  is  unbroken,  comprising  cloud  of 
variable  optical  depth  whilst  the  bottom  layer  is  broken 
with  cloud  of  uniform  optical  depth  (after  Platt,  1983). 


> 

t 

Table  6.  Parameters  emp loyea  in  sins le  c loua  layer  ca leu lations.  T_c 

and  T„  are  the  cloud  and  surface  temperatures  respectively  g  is  the 
£ 

surface  a  1  bedo,  represents  the  incidence  ang le  of  solar  radiation  at 
the  cloud  whilst  §  denotes  the  nadir  angle  o f  the  cloud  at  the 
sate  1 lite*  Q)  is  the  cloud  albedo  as  a  function  of  its  visible 

optical  depth  and  so..l  a_r  radiation  incidence  ang  le  whilst  &  i_s  the 
ratio  of  visible  optical  depth  divided  by  the  infrared  absorption 
optical  depth  (after  Platt,  1  9b3  )• 


Boundary  layer  cloud 


high  cloud 


each  other*  In  the  latter  case  the  shape  of  the  curves  is  determined 
by  four  independent  variables  namely,  the  cloud  fraction  ana 
emissivity  of  the  two  overlapping  layers.  Figures  25(b)  and  25(c) 
illustrate  curves  depicting  typical  atmospheric  cloud  systems.  The 
basic  equations  have  been  simplified  here  for  cases  of  cirrocumulus 
overlying  altocumulus  (Figure  25(b))  and  cirrostratus  overlying 
cumulus  (or  altocumulus)  (Figure  25(c))*  Again  variations  in  layer 
cloua  fraction  and  cloud  emissivity  are  reflected  in  albedo  ana 
brightness  temperature  changes.  In  addition  Platt  (1983)  illustrates 
the  possible  effects  due  to  infrared  scattering  when  the  cloua  layers 
become  deeper  ana  the  model  assumption  of  an  isothermal  cloud  becomes 
invalid. 

The  model  curves  have  been  tested  against  multiple  cloud  layer 
systems  over  ocean  and  land  surfaces.  Albedo/brightness  temperature 
histograms  are  constructed  from  satellite  measurements  and  model 
curves  fitted  to  assess  the  level  of  agreement.  Figure  26  shows  the 
result  of  fitting  a  model  curve  for  semitransparent  cloua  to  a 
bispectra  1  histogram  for  data  from  a  complex  tropical  convective 
system.  The  consistency  is  evident  along  much  of  the  histogram's 
major  axis  auguring  well  for  the  model  structure  in  this  case. 

The  overall  value  of  this  radiative  transfer  algorithm  is  so  far 
unknown.  More  comprehensive  model  tests  are  still  necessary, 
particularly  on  convective  and  frontal  clouds,  as  well  as  on  low  level 
stratus  which  require  the  inclusion  of  features  such  as  anisotropic 
cloua  reflection  into  the  model.  Platt  (19&3)  reports  that  any  future 
development  leading  to  possible  incorporation  into  the  ISCCP  algorithm 
will  have  to  follow  identification  of  a  successful  method  for 
automating  the  removal  of  cloud  parameters  from  the  input  bispectral 
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Figure  26 
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Two-dimensional  bispectral  histogram  for  a  500  km  x  500 
km  area  containing  tropical  clouds  on  11  April  I960, 
latitude  2.5°N,  longitude  137.5°E»  The  lines  represent 
two  predicted  bispectral  curves  for  a  single  layer  of 
semitransparent  cloud  with  cloud  top  temperatures  of  200 
and  210  K,  respectively  (after  Platt,  1985)* 
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3  Feasibility  of  Using  Surface  Observations  < 

Throughout  the  discussion  of  retrieval  algorithms  in  Section  2  j 

the  need  for  improved  development  ana  subsequent  testing  has  been  ! 

i 

i 

emphasised  because  many  of  the  methods  are  still  in  their  infancy  and  j 

have  been  tested  only  on  limited  regions  on  a  limited  time  scale  in 
the  course  of  which  the  full  range  of  conditions,  heterogeneous 
surface  properties,  time-varying  surface  problems  and  all  the  major 
global  cloud  systems  have  not  been  encountered.  The  aim  of  those 
involved  in  the  1SCCP  pilot  study  in  1 9 6 1  was  to  try  to  rectify  this 
in  part  but  the  caveats  still  apply  to  developments  that  have  taken 
place  since  (spatial  coherence,  maximal  clustering,  bispectral  ] 

curves).  Perhaps  the  most  pertinent  problem  to  climate  modellers  at 
present  is  the  lack  of  a  uniform  global  cloud  climatology.  The  data 
set  that  will  hopefully  emerge  from  ISCCP  will  be  wholly  satellite 
based.  The  retrieval  algorithm  that  was  proposed,  after  the  pilot 
study  was  not  specifically  any  of  those  tested  but  combined  the  'best' 

features  of  some  of  the  individual  methodologies.  There  will, 

however,  still  remain  the  problem  of  how  to  validate  the  data  set. 

The  most  readily  available  method  (ana  perhaps  the  only  conceivable 
one)  lies  in  comparison  with  surface  synoptic  reports  but  the  paucity 
of  such  exercises  (Reynolds  ana  Vonaer  haar,  1977)  had  until  recently 
left  something  of  a  void  in  this  area. 

The  problems  of  satellite  retrieval  are  now  fairly  well 
established  (Rossow  et  a 1«,  1965)  and  the  question  arises  can  surface 

observations  be  added  to  or  in  some  way  calibrated  to  satellite 

retrievals  so  as  to  improve  the  quality  of  the  latter.  The  problem  is 
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addressed  by  Sfeze  ejt  ajU  (1967)  which  describes  details  of  a  20  oay 
period  m  the  summer  of  1983  during  which  satellite  to  surface  cloud 
amount  comparisons  were  made  over  Western  Europe  using  data  from 
surface  stations  in  Britain  and  France  and  the  corresponding  hETEGSAT 
imagery.  The  exercise  formed  a  part  of  the  1SCCP  preliminary 
validation  exercise  and  was  colocated  with  one  of  the  special  study 
regions  designated  by  ISCCP.  hethoos  of  the  surface  retrieval  are 
given  in  Sfeze  e_t  a  1 .  (  1  987)  whilst  the  satellite  algorithm  was  the 
four-parameter  (Ih/IRva„/V IS/ V ISva„)  clustering  technique  (Desbois  and 
Sfeze,  1984a). 

Overall,  the  level  of  agreement  in  total  cloud  amount  between 
satellite  and  surface  observer  was  +  1  okta  (or  +12%)  for  64%  of  all 
cases.  The  study  raises  many  interesting  and  pertinent  points  which 
are  considered  briefly  here. 

beginning  with  total  cloud  amount,  there  was  found  to  be  greatest 
agreement  in  either  clear  sky  situations  or  single-layered  complete 
overcasts  but  low  cloud  amounts  (important  for  considering  surface 
radiative  fluxes)  were  considerably  underestimated  by  the  satellite 
whenever  middle  or  upper  level  cloud  overlaid  the  low  cloud  deck. 
This  is  possibly  a  case  in  which  surface  reports  could  assist  in 
correcting  satellite  measurements  (layer  obscuration). 
Correspondingly  it  was  found  that  the  surface  observer  underestimated 
high  cloud  amount  (>6  km  altitude)  when  it  was  obscured  from  view  by 
lower  decks.  In  this  situation  the  success  with  which  the  algorithm 
was  found  to  detect  cirriform  clouds  well,  becomes  important.  In  this 
study  serious  discrepancies  arose  for  such  cases,  which  were  directly 
related  to  the  success  with  which  the  algorithm  could  detect  cirrus 


a  wide  range  of  alternative  algorithms. 

The  use  of  METEGSAT  radiances  at  middle  latitudes  raises  problems 
of  viewing  geometry  and  the  pixel  resolution  of  6  km  x  fa  km  means  that 
many  pixels  may  be  partially  cloud  covered.  At  the  latitudes  concerned 
M  ETEOSAT  will  view  the  sides  of  clouds,  particularly  those  of 
developing  cumuliform  clouds  which  were  observed  to  occur  ana  small 
gaps  in  layers  will  also  be  missed.  however,  perspective  problems  are 
not  confined  to  the  satellite  alone.  The  surface  observer  also  fails 
to  detect  gaps  in  the  cloud  cover  at  low  elevation  angles  (Malberg, 
1973)  and  tends  to  overestimate  cloud  amount  near  to  the  horizon 
(Merritt,  1966).  The  perspective  effects  can  be  shown  to  be  of  similar 
magnitude  for  both  surface  and  satellite  retrievals  but  should  not  be 
disregarded  altogether. 

The  representativeness  of  the  surface  stations  is  another 
important  area.  In  an  ideal  comparison  a  spatially  homogeneous,  cense 
observing  network  could  be  colocated  with  the  satellite  image  pixels. 
Reality  is  far  removed  however,  and  considerations  to  the  area 
representative  of  each  surface  station  are  given  in  Stee  eX  a  1.  (196b, 

1 9&7 )•  It  should  also  be  conceded  that  a  surface  observer  is  an 
essentially  subjective  analysis  instrument.  Perhaps  the  use  of  all¬ 
sky  cameras  could  provide  a  more  objective  retrieval  (Ackerman  ano 
Cox ,  1961). 

The  question  of  incorporating  surface  observations  into  the 
satellite  retrieval  remains  a  very  difficult  one,  A  method  of 
combining  the  two  data  types  realistically  has  yet  to  be  proposed. 
One  possibility  concerns  the  work  of  Hahn  et  a  I.  (19691  which  consists 
of  surface-derived  cloud  climatologies  expressing  the  probability  that 
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given  the  presence  of  one  cloud  type,  another  type  is  present.  If  the 
satellite-derived  cloud  type  is  assumed  to  be  the  referenced  type  the 
extra  information  provided  by  these  climatologies  could  improve  future 
retrievals. 

Much  work  remains  to  be  done  on  algorithm  testing,  improvement 
and  validation.  It  now  appears  certain  that  any  future  global 
retrieval  algorithm  will  be  a  heterogeneous  composite  of  a  range  of 
satellite-based  techniques.  It  is  therefore  reasonable  to  propose  the 
inclusion  of  surface  observations  where  they  would  improve  the  cloud 
characterisation. 

Summary  and  Conclusions 

This  report  describes  the  first  phase  of  a  three  year  study  of 

cloud  retrieval  techniques.  The  aim  this  year  has  been  to  review 

existing  techniques  with  particular  attention  being  paid  to  satellite- 

basea  retrievals  and  to  novel  proposals  made  since  the  ISCCP 
h 

intercomparison  was  initiated  in  1981/82.  The  somewhat  unsatisfactory 
conclusion  that  must  be  drawn  is  that  while  many  new  techniques  offer 
useful,  even  uniquely  successful,  retrieval  capabilities  in  certain 
well-defined  situations  they  have  not  yet,  and  perhaps  can  never,  be 
shown  co  be  applicable  to  other  areas,  climates  ana  configurations. 

Two  recommendations  are  made:  (i)  there  should  be  a  broader 
"validation"  of  many  of  the  techniques  reviewed  here  and  (ii)  an 
attempt  should  be  made  to  incorporate  surface-basea  observations  in 
retrievals.  It  seems  likely  that  in  many  situations  surface-basea 
observations  offer  the  only  satisfactory  means  of  detecting  (a)  small- 
scale  structure,  especially  fair  weather  cumulus  and  gaps  in  cloud 
decks,  (b)  thin  cirrus  and  (c)  low  cloud  amount  and  type  in  the 
presence  of  miaale  ana  high  cloua  layers.  These  preliminary 


conclusions  seem  to  justify  the  U.S.  Air  Force's  use  of  "all  available 
data"  in  their  construction  of  global  nephanalyses. 

Presentations  and  Publications 

As  most  of  this  year's  work  has  been  of  an  introductory  nature 
consisting  primarily  of  literature  review  no  publications  have  ensued. 
Two  presentations  have  been  made: 

Postgraduate  seminar  series  at  Liverpool  University  on  15th  June  1 9b 6 
entitled 

"Automated  analyses  of  satellite  radiance  data". 

Seminar  at  NATO  School  on  Remote  S_e  n  s^  _i  n  £  i^n  Meteorology  and 
Climatology  in  Dunaee  (Aug.  /  Sept.  1966)  entitled 
"A  brief  introauction  to  satellite  cloud  retrieval" 
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